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Abstract. We present the combined results on electron-pair production in 158GeV/n Pb-Au 
{\/^~ 17.2 GeV) collisions taken at the CERN SPS in 1995 and 1996, and give a detailed account of 
the data analysis. The enhancement over the reference of neutral meson decays amounts to a factor of 
2.31±0.19 {stat.) ± 0.55 (syst.) ± 0.69 (decays) for semi-central collisions (28% a/ageo) when yields are in- 
tegrated over m > 200MeV/c^ in invariant mass. The measured yield, its stronger-than-linear scaling with 
Nch, and the dominance of low pair pt strongly suggest an interpretation as thermal radiation from pion 
annihilation in the hadronic fireball. The shape of the excess centring at m ~ 500 MeV/c^, however, cannot 
, ^ , be described without strong medium modifications of the p meson. The results are put into perspective 

' by comparison to predictions from Brown-Rho scaling governed by chiral symmetry restoration, and from 

the spectral- function many-body treatment in which the approach to the phase boundary is less explicit. 
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1 Introduction 

Strongly interacting matter under extreme conditions of 
temperature and density is being created by colliding heavy 
nuclei in fixed target experiments at the Super-Proton 
Synchrotron (SPS) at CERN and at the Relativistic Heavy- 
Ion Collider (RHIC) at BNL. The motivation derives from 
the quest to discover 'quark matter', the quark gluon plas- 
ma (QGP), in which quarks and gluons are deconfined 
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and chiral symmetry is restored 0, two complementary 
facets of the phase transition predicted by Quantum Chro- 
modynamics (QCD). This transition is expected to occur 
at a critical energy density of e «0.7 GeV/fm^ and a 
temperature of 170-180 MeV as finite-temperature lattice 
calculations have shown The initial energy density at 
full SPS energy is appreciably larger and reaches about 
3 GeV/fm^ in central Pb-Pb collisions, adopting Bjorken's 
longitudinal expansion scenario for a simple estimate 
While the existence of the QGP has not yet been proven, 
there is circumstantial evidence that its transient forma- 
tion is imprinted, in one way or another and to varying 
degrees, on measured distributions of final-state hadrons. 
Since the fireball terminates in an exploding multi-hadron 
final state, QGP signatures, if not collective in character, 
are prone to be masked by hadronic interactions. 

Photons and dileptons are potentially more direct pro- 
bes of the early collision stages since they escape from 
the impact zone nearly undisturbed by final-state inter- 
actions and have their largest emission rates in hot and 
dense matter. Moreover, according to the vector domi- 
nance model 0], dilepton production is mediated in the 
hadronic phase by the light neutral vector mesons p, w, 
and (j) which mark the low-mass region by distinctive res- 
onance peaks. Among these, especially the short-lived 
(0(770) meson (t = 1.3 fm/c) has acquired a key role as test 
particle for 'in-medium modifications' of hadron proper- 
ties close to the QCD phase boundary |^||^||7). Changes in 
position and width of the p have been advocated already 
20 years ago as precursor signals of the chiral transition 8 . 
Restoration of chiral symmetry in hot and dense matter 
has become one of the heavily discussed and exciting is- 
sues in non-perturbative QCD thermodynamics |n]l7|l9|llOI 
ITT] as the melting of the chiral condensate should cause 
rather drastic changes of the properties of the light vec- 
tor mesons and thereby on the structure of dilepton spec- 
tra. Spectral function calculations on the lattice seem still 
far from providing model-independent guidelines for the 
study of thermal modifications of hadron properties . 

An enhancement of low-mass lepton pair production 
was first reported by CERES [TWr^ and HELIOS-3 [H 
Uni on the basis of 200 GeV/n S-Au and S-W data, re- 
spectively. The CERES e+e~ mass spectrum is shown in 
Fig. 11.11 together with the reference spectra for 450 GeV 
p-Be and p-Au collisions While the p-A data are 

reproduced within errors by final state Dalitz and direct 
decays of neutral mesons as known from p-p collisions, 
electron pairs from S-Au collisions reveal a substantial en- 
hancement in the mass range above 250 MeV/c^. 

At top SPS energy and close to the critical tempera- 
ture, the prime candidate for 'thermal radiation' from the 
hadronic phase of the fireball |17II18| is pion annihilation, 

vr^TT^ ^ p — > e^e^. (1-1) 



^ The notation of the ordinate in Fig. II. H is no longer in use; 
it should read as in all other mass spectra shown in this paper. 
The 'cocktail' has received minor adjustments in the meantime 
which do not affect any of the conclusions drawn. 



This thermal process with a threshold at 2 m^r is dynami- 
cally enhanced via the electro-magnetic form factor of the 
pion by the p resonance ,19, , with a dilepton branching of 
1 in 10'^. Yet, the p serves not only as test particle, but its 
strong coupling to the tttt channel makes it also a major 
constituent of hot hadronic matter. Numerous theoretical 
approaches incorporating pion annihilation using vacuum 
properties of the p meson, failed without exception to de- 
scribe the data PU] . 

This suggested in-medium changes of the p spectral 
function that shift dilepton strength down to lower masses. 
The first calculations that were successful in describing 
both the CERES and the HELIOS-3 data made use of 
the scaling conjecture of Brown and Rho 5^, which postu- 
lates that the mass of non-strange vector mesons decreases 
in dense matter together with the scalar quark conden- 
sate, the order parameter of the chiral transition. This 
'dropping mass' scenario received independent support by 
work on QCD sum rules (21] . and it turned out like a 
tailor-made concept: linked to a fireball model |22ll23j or 
embedded into transport calculations |24II25| . it gave ex- 
cellent fits to the data. 

But it was also pointed out |26II27II28II29| that chiral 
symmetry considerations alone would only require that 
masses of chiral partners, here the vector meson p and 
the axial vector meson oi, become degenerate, but by no 
means necessarily massless. Moreover, significant mixing 
must accompany any mass shift when approaching the 
phase transition along T 317. The 'dropping mass' scal- 
ing idea has been very recently revisited by the original 
authors |Zj welcoming an alternative scheme of how chiral 
symmetry might be restored, the Georgi vector limit j31| 
in which the chiral partner of the (longitudinal) p is the 
Goldstone pion, both becoming massless in approach of 
the chiral transition |32II33| . Thermal modifications of hadron 
properties in general, and of in-medium spectral functions 
in particular, have not yet come within reach of QCD lat- 
tice calculations with finite baryon density 2 . Very pre- 
liminary results suggesting dropping vector meson masses 
have been reported 

An alternative approach to explain the low-mass dilep- 
ton enhancement in the CERN SPS data focused on the 
calculation of spectral functions in a hot and strongly in- 
teracting hadron resonance gas at finite baryon density 
by conventional many-body techniques f5ll35j . These con- 
firmed earlier calculations of the two-pion self-energy in 
nuclear matter which indicated that the p spectral func- 
tion suffers significant broadening but only negligible shift 
in mass 36 , 37 . The many-body calculations combined 
with suitable reaction models describe the observations 
actually quite well |38II39II40II41| : the low-mass wing of the 
broadened in-medium p spectral function receives strong 
thermal Bose enhancement which results in an amplified 
dilepton strength considerably below the vacuum p posi- 
tion, while simultaneously the yield at the vacuum po- 
sition was depleted, consistent with approximate unitar- 
ity [ig. 

Our discussion so far was limited to aspects of dilepton 
emission rates with possible modifications by the medium. 
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Fig. 1.1. CERES inclusive e+e" mass spectra of 450 GeV p-Be, p-Au, and 200 GeV/n S-Au collisions |16I12| . Plotted is the 
number of electron pairs per charged particle, both in the acceptance and per event. Contributions from various hadron decays 
as expected from p-p collisions are shown together with their sum (thick line), and the systematic error on the latter is indicated 
by the shaded area. In S-Au a higher single-electron pt cut reduces the 7r° Dalitz component. See text and footnote. 



However, total pair yields derive from space-time integra- 
tion over a priori unknown density and temperature pro- 
files which are usually modelled by hydro-dynamical |48l 
or microscopic transport 35,23,24 calculations, or 
fireball models 22 ,45 . Certainly, the external inputs, e.g. 
to the hydro-dynamical and fireball calculations, have to 
conform with whatever knowledge there is on initial condi- 
tions, depending on collision geometry, and on the {fj,, T) 
coordinates of the trajectory in the phase diagram. The 
medium modifications of the p suggested by the SPS re- 
sults seem to require a strongly interacting, hot and dense 
hadronic fireball with sufficient time spent between hadro- 
nisation and thermal freeze-out. If this time would be 
insufficient, the enhanced dilepton production may have 
stronger links to the hadronisation stage or the plasma 
phase than hitherto assumed. 

The CERES Collaboration measured dilepton produc- 
tion in 158GeV/n Pb-Au collisions in 1995 and 1996 with 
a greatly improved setup 46 , 47 compared to the sulfur- 
beam experiment. The main objective of the Pb runs has 
been achieved: to corroborate with a large statistics sam- 
ple the enhanced dielectron production at low masses for 
the heavy Pb-Au collision system 47,48,49,50,51 . Im- 
proved background rejection was achieved which was com- 
pulsory in an environment of very large rapidity density 
of hadrons and secondary photons. Among the physics 
goals considered most important for further insight into 
the nature of the processes at work was the centrality de- 
pendence of the enhancement. 

The significance of the baryon chemical potential for 
in-medium modifications at SPS energies prevailing over 
that of pion number, or temperature, as found in most cal- 
culations j35ll()] , but not in all "^2] , was given experimental 



support by the recent finding of an even somewhat larger 
enhancement measured in the CERES Pb-Au low energy 
run at 40 GeV/n [SSIISHI, compared to that at 158 GeV/n. 

This paper presents the combined results of all data 
on electron pair production in 158 GeV/n Pb-Au collisions 
taken by the CERES Collaboration in the years 1995 and 
1996. Most of the analyses were performed in the course of 
Doctoral Dissertations in Heidelberg, Darmstadt and Re- 
hovot dealing with the 1995 and 1996 |KKII^?)IIK7| 

data. Publications of analysis results of the 1995 data 03 
48,49 and the 1996 data 50,5^ are superseded by the 
combined results presented here. There are no major de- 
viations of the unified results reported in this paper to 
those pubfished previously for the separate data sets. 

The paper starts with a description of the experimen- 
tal setup and summarises the instrumental means CERES 
has at its disposal to cope with background. A detailed 
description of the data analysis is given in sect. |31 which 
concludes with the centrality determination. The Monte- 
Carlo simulation method applied for measuring reconstruc- 
tion efficiency and optimising the rejection of combinato- 
rial background is addressed in sect. 01 The 'cocktail' of 
hadron decays which serves as an important reference for 
electron-pair production is discussed in sect.[S| Results of 
both data sets are presented in sect. El which also includes 
a discussion of statistical and systematic errors. The final 
mass and transverse momentum spectra are presented and 
compared to the hadronic cocktail in sect. [3 Section|Hlcon- 
tains a physics discussion, from an experimentalist point 
of view, on the comparison of data to current theoretical 
models. The paper concludes by summarising what has 
been achieved and which issues are still open but might 
be clarified in the not too distant future. 




Fig. 1.2. CERES setup during the 1995/1996 Pb-beam runs 



2 The CERES experiment in 1995/96 

CERES is dedicated 59 to the measurement of electron 
pairs in the low-mass range from 50 MeV/c^ up to about 
~1.5 GeV/c^; an upper mass limit is imposed by count- 
ing statistics due to the rapid decline in cross section. 
The spectrometer covers the pseudo-rapidity region close 
to mid-rapidity, 2.1 < < 2.65. It is axially symmetric 
around the beam and has 27r azimuthal coverage. Trans- 
verse pair momenta are accepted down to ~20 MeV/c for 
masses above ~400 MeV/c^. These are great assets when 
investigating soft processes. CERES maintains transpa- 
rency for hadrons and photons as strictly as possible. 

A schematic view of the spectrometer in the run peri- 
ods 1995 and 1996 is shown in Fig. II. 21 At the heart of the 
spectrometer are the two coaxial ring-imaging Cherenkov 
detectors HOj, one (RICH-1) within the other (RICH-2) 
along the beam and separated by a compact super-cond- 
ucting solenoid for momentum analysis. A doublet of sili- 
con-drift chambers (SiDC) [HI] replaces the previously 
used single SiDC to enable precise charged-particle track- 
ing into RICH-1 for rejection of close tracks, and for off- 
line measurement of charged multiplicity. A multi-wire 
proportional counter (Pad Chamber) with pad readout 
was added behind the mirror of RICII-2 to provide exter- 
nal tracking downstream of RICII-2. The new tracking de- 
tectors outside the field were added to cope with the high 
multiplicities of Pb- Au collisions I62ll63j . These upgrades 
are described in Ref. fl6 and Ref. [HI] for RICH and 
SiDC, respectively. A multiplicity detector (MD) of plastic 
scintillators behind the Pad Chamber serves as first-level 
trigger device. Below we introduce the individual detector 
components in order of their arrangement along the beam. 



2.1 The target area 

The target area, hardly visible in Fig. 11.21 and enlarged 
in Fig. 12.11 comprises the segmented target, the doublet 



of SiDC's, the light-collecting parts of the interaction- 
vetoing beam counter BC3 (see sect. I2.6|l and interfaces 
to the adjacent parts of the spectrometer. It is housed 
within a hollow cylindrical recess into RICH-1 which is 
lined by a cylindrical tungsten mantle of 20 mm thick- 
ness which shields the UV-detectors from heavily ionis- 
ing particles emerging in backward direction from the tar- 
get. The recess is hermetically closed towards the radiator 
of RICH-1 by a double- walled aluminised mylar window 
(2 X 50 /im), the intermediate volume being ventilated 
with nitrogen. The SiDC's and a thin mirror for BC3 are 
mounted within a double- walled tube of 160 mm diameter, 
made of aluminium and carbon fibre which is water cooled 
and provides a laminar flow of dry nitrogen for cooling 
the SiDC's and the front-end chips mounted on the same 
mother board. A light guide transfers the Cherenkov light 
received from a thin mirror. 



UV ]ight-<juxie cabfe 



Beam 










1 1 1 1 r M r 




Vacuum 


Au target 






Fig. 2.1. Schematic cross section of the target area along the 
beam direction showing the segmented target, the doublet of 
SiDC's, and in between the small mirror of BC3 which directs 
Cherenkov light from every passing Pb ion into a light guide 
for vetoing the interaction trigger. The acceptance of the RICH 
detectors is indicated (dashed lines). 
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The beam enters the target area via an evacuated Al 
tube that reaches until a few millimetres short of the seg- 
mented target and is sealed by a thin mylar window . 

The segmented Au target consists of 8 discs of 25^m 
thickness and 600/Ltm diameter spaced uniformly by 3.2 mm 
along the beam. The total target thickness of 200/im Au 
corresponds to an interaction length of A/ A/ = 0.83% for 
Pb-Au collisions. The segmentation assures that only a 
fraction of the total radiation length is effectively seen by 
photons and electrons propagating into the spectrometer 
acceptance. Under worst conditions of full illumination, 
there is only a chance of about 20% that the next down- 
stream disc is within the acceptance. The effective radi- 
ation length is X/Xo ~ 0.55%, compared to 3% for an 
unsegmented target of identical interaction length. 

The Au discs are supported by 2.5 /im thick mylar 
foils, and the entire assembly is contained in a thin-walled 
(0.5 mm) carbon fibre tube of 60 mm diameter; the Au 
discs are accurately aligned to the axis of the tube. The 
Au discs contain 90% of the beam over the full length 
of the target once collimators and magnet settings of the 
beam line have been optimised. The beam diameter was 
measured to have a Gaussian envelope of cr= 220 /im. 

2.2 The silicon-drift telescope 

The SiDC's fulfil several purposes in the overall concept: 
(i) to locate the interaction vertex within the segmented 
target, (ii) to maintain a precise two-point charged-particle 
tracking free of ambiguities, improving thereby the mo- 
mentum resolution of the spectrometer, (iii) to assist in 
the ring pattern recognition of RICH-1, and (iv) to reject 
conversions and Dalitz pairs which are not resolved in the 
RICHes, by detecting a double pulse-height signal or two 
close tracks. 

The SiDC telescope made up of two closely spaced 
cylindrical SiDC's was the major element in upgrading 
CERES for the Pb-beam experiments 63. CERES is the 
first experiment which successfully implemented 65 this 
detector concept 66^, and it did so from the start, how- 
ever, with only a single detector j61ll67j . This allowed to 
check that the electron pair originates from a common 
vertex, but was of course insufficient for full tracking. 

Since cylindrical SiDC's seem not very well known, we 
insert here a short description. The active area is prac- 
tically the full area of a 3-inch-diameter wafer 280 /im 
thick which has a central hole of about 6 mm diameter 
for the passage of the beam. The signal electrons gener- 
ated by ionisation of charged particles traversing the wafer 
drift radially outward to an array of 360 anodes located 
at the periphery of the detector. The radial coordinate r 
of the point where a charged particle crossed the detec- 
tor plane is measured by the drift time t of the electron 
cloud. The charge sharing between neighbouring anodes 
measures the azimuthal coordinate cf). The pair of coor- 
dinates (r, (jj) is provided for each crossing charged parti- 
cle for events with a total charged multiplicity of several 
hundred. The longest drift distance is about 3 cm. The 
nominal value of the drift field of 500 V/cm results in 



a maximum drift time of about 4 /is. The drift field is 
provided by means of 240 concentric electrodes with 
130 fim pitch on both sides of the detector, suitably bi- 
ased by an implanted voltage divider. A major innovation 
in the design was a 'sink anode' providing a path for the 
leakage current generated at the Si-Si02 interface away 
from the signal anode, i.e. without contributing to the an- 
ode leakage current. 

The ideal circular electrode shape could not be de- 
signed in early 1990 due to limitations in the software 
controlling mask production. Electrodes were shaped in- 
stead as regular polygons of 120 sides. We were surprised 
to see that the distribution of hits displayed peaks at every 
third anode 68 . The 'efficient' anodes were those located 
in the central part of each 3° triangle forming the polygon. 

The telescope implemented in 1995 consisted of a 3- 
inch detector followed by a 4-inch detector, the latter of 
the novel AZTEC design which eliminated the focus- 
ing problem. The spacing of the two detectors was 14.3 
(15.0) mm and the distance from the target centre to 
SiDC-1 amounted to 98.5 (110) mm in 1995 (1996). 

In production of the 3-inch detector^, the front and 
back lithographies were rotated with respect to each other 
by 1.5° to reduce the non-radial components in the drift 
field. This trick reduced to negligible levels the focusing 
effect of the central anodes that had been a severe obstacle 
for reaching the design azimuthal hit resolution |68j . 

In 1996, both detectors were of the 4- inch type^. The 
larger anode radius of 42 mm allowed to increase the dis- 
tance to the target for lower hit occupancy, especially at 
small radii. The sensitive area of the detectors is increased 
from 32 cm'^ for the 3-inch devices to 55 cm^ for the 4- 
inch detectors. The 'field cage' polygons consist of 277 
concentric p"*" electrodes each having 360 instead of 120 
sides before, and the respective p+-implantations on the 
two wafer sides are rotated by 0.5° with respect to each 
other for near perfect radial field geometry. Another novel 
property of the AZTEC detectors is the interlaced anode 
structure: each one of the 360 anodes is subdivided into 5 
segments, two of which (extending over 16% of the anode 
pitch of 1 degree) are interlaced to the closest of the neigh- 
bouring anodes to enforce charge sharing. This provides 
a more accurate azimuthal position measurement when 
calculating the centre of gravity of the distribution. Com- 
pared to the '95 runtime, the new design improves the 
azimuthal resolution from 2.5 to 1 mrad. The resolution 
in radial direction is 30 /im both for the 95 and 96 set-ups. 

Charge signals from 3-inch detectors are amplified by 
32-channel front-end OLA chips placed on the detector 
motherboards which were developed to test ALICE silicon- 
drift prototype detectors [^21 and were produced in a cus- 
tom bipolar process'*. They consist of a charge-sensitive 



^ produced by SINTEF, 0134 Oslo 3, Norway 
^ produced by EURISYS Mesures, F-67380 Lingolsheim, 
France 

* Owned by Tektronix at the time of production; since then 
Maxim Integrated Products 
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preamplifier, a quasi-Gaussian shaper^ and a symmetrical 
line driver. The rather high gain (30 mV/fC) gave rise to 
wild collective oscillations. Stable operation was achieved 
only with the introduction of 20 f2 damping resistors con- 
nected in series into the 50 12-terminated output lines. 

For the readout of the 4- inch detectors, new 16-channel 
front-end chips had been designed along a CMOS con- 
cept |7()| which incorporated bipolar drivers and were well 
adjusted in shaping time (37 ns), gain (9 mV/fC), dy- 
namic range (5 mips^), and low equivalent noise charge 
(140 e~) to meet our requirements.^ 

An outside buffer stage transmits the bipolar signals 
over 40 m flat cables to the FADC's {flash analog to digital 
converters)^ which sample the data with 50 MHz and store 
it 256 bytes deep. This corresponds to a drift time range 
of 5.12 fj,s. Digitisation is 6 bit with non-linear character- 
istics.® Data are continuously sampled until interrupted 
by an external 'stop' signal so that the channel memory 
contains always the last 5.12 /is of data. 

Since the trigger signal and clock are asynchronous, 
there is a random phase difference producing a time jitter 
of 20 ns/-\/T2 rms. It is measured with a TDC {time-to- 
digital converter) , and correction is done off-line. 

Each of the four FADC crates per detector houses a 
SIM {scanner interface module) which scans the data af- 
ter the trigger was received for contents above a prede- 
fined threshold {readout threshold). Readout is activated 
whenever the threshold is surpassed in two successive time 
bins, and stopped, if contents in two successive time bins 
fall below it. The contents of five preceding channels are 
also readout {pre-samples) for off-line reconstruction of 
the baseline. 



2.3 The RICH detectors 

The radiators are filled with methane at atmospheric pres- 
sure. The high Cherenkov threshold of jthr — 32 en- 
sures that more than 95% of all charged particles pass 
without creating Cherenkov light ('hadron blind track- 
ing'). The Cherenkov light is reflected backward onto 2- 
dimensionally position-sensitive gas detectors which are 
separated from the radiator volume by UV-transparent 
windows. By their upstream position with respect to the 
target, the UV-detectors are not exposed to the huge for- 
ward flux of charged particles. The price to be payed for 
this geometry is the limited acceptance in polar angle 0, 
indicated by the lines in the upper part of Fig ll.2l 

High-energy electrons produce Cherenkov rings with 
asymptotic radius, Roo= ^/ithr — 30 mrad. The differ- 

^ The time constant t— 38 ns, about twice the design value, 
deteriorated the potential double-pulse resolution, but avoided 
a large ballistic deficit over the full drift. 

^ minimum ionising particles 

^ produced by AMS in 0.8^m biCMOS technology. 
® Series DL300 of Fa. B. Struck, Tangstedt near Hamburg 
3 Channel.No.(0...63) = 256 ■ Uin/(0.2 -\- 3 • Uin), input volt- 
age Uin in Volt. 

^° different in RICH-1 and RICH-2 by about one unit. 



ence in radiator lengths (86 cm and 175 cm for RICH- 
1 and RICH-2, respectively) is partially compensated by 
better UV transmission in RlCH-1 (CaF2 window) com- 
pared to RICH-2 (quartz window), so that the asymptotic 
number of photons per ring, 10.8 and 11.5, for RICH-1 and 
RICH-2, respectively, come out rather similar^^. By the 
same reason, photon detection in RICH-1 reaches farther 
into the UV. 

In both RICHes spherical mirrors focus the Cherenkov 
photons radiated from a straight trajectory back onto a 
ring image in the focal plane of the UV detectors. As the 
mirror in RICH-1 is traversed by all electrons before the 
second ring image for momentum measurement is taken 
in RICH-2, there are stringent physics reasons to keep the 
radiation length as low as possible: besides reducing the 
number of external conversions, it is the multiple scatter- 
ing of low-momentum electrons which reduces the detec- 
tion efflciency for soft pairs and deteriorates the momen- 
tum resolution. The mirror of RICH-1 therefore is made 
very thin (1.1 mm) so that it adds only 0.4% of a radia- 
tion length. It is based on a laminated carbon fibre struc- 
ture which defines the spherical geometry^^. An evapo- 
rated coating of aluminium protected by magnesium fluo- 
ride achieved persistent UV reflectivity of 80% at 300 nm. 

The UV-detectors consist of a conversion space fol- 
lowed by two parallel-plate avalanche stages and a multi- 
wire proportional chamber. The originally planned mode 
of running only with parallel-plate amplification was aban- 
doned in favour of an added multi-wire stage, following a 
painful learning process on spark break down problems in 
pure parallel-plate schemes 71 . 

The operating gas is He -I- 6% CH4 at atmospheric 
pressure -I- TMAE-saturated^'^ vapour at 40°C as photon 
converter. The use of TMAE demands that the UV detec- 
tors be kept hot to avoid condensation. To avoid temper- 
ature gradients across the delicate UV-transparent win- 
dows separating the detectors from the radiators, the en- 
tire spectrometer is kept hot at about 50° C. The UV 
detectors operate at a total gain of about 2-10^ for high 
photon detection efficiency (~ 85%). The ion clouds pro- 
duced in the last wire amplification stage induce signals 
on the pads. The latter form a grid of pitch 2.74 mm and 
7.62 mm, and the resulting total number of pads is 53,800 
and 48,400 in RICH-1 and RICH-2, respectively (72]. 

The UV detectors have been operated without open- 
ing since 1991. During the 1995 run, UV-1 degraded in 
performance. The detector could not be operated at the 
desired gain of 2x10^ without an excessive spark rate. 
Early in 1996, the UV-1 detector was opened and all mesh 
electrodes, in particular the cathodes and the multi-wire 
plane showed some kind of deposit. Most mesh electrodes 
were exchanged and the wire anode subjected to ultra- 
sonic cleaning. The refurbished detector performed very 
well during the 1996 run |55j . 

The numbers of Ref. |60l measured with ethane are in- 
creased by one photo-electron due to the larger bandwidth in 
methane. 

manufactured by MAN Technologic AG 
" Tetrakis-di- Methyl- Amino- Ethylen 
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2.4 Deflection in the magnetic field 

The magnetic field for momentum analysis is generated by 
two super-conducting solenoids carrying currents in oppo- 
site sense. Charged particles experience an azimuthal de- 
flection between the two RICHes which is inversely pro- 
portional to the momentum, 

and the sense of which, for fixed polarity of the field, de- 
fines the charge sign. The constant is 4>o= 146 mrad GeV/ c. 
Between SiDC and the Pad Chamber the deflection is only 
66% of this value, (f)o— 96 mradCeV/c.^'* To first order, 
the polar angle is not affected. Particles deflected in az- 
imuth by Z\0 encounter a small second-order deflection 
towards the beam axis which amounts to 

Ae = -78 [A^^f . (2.2) 

Two sets of warm correction coils are tuned to achieve a 
field- free radiator in RICH-1 and to align the field lines in 
the radiator of RICH-2 parallel to the particle trajectories 
from the target. This way straight trajectories inside both 
radiators are achieved. Moreover, the absence of deflection 
in the flrst RICH detector allows to identify conversion 
and Dalitz pairs by their small opening angles. 

2.5 The Pad Chamber 

The last tracking detector is a multi-wire proportional 
chamber located closely behind the mirror of RICH-2, at 
a distance of about 3.3 m from the target. It has an in- 
ner and outer radius of 42 cm and 85 cm, respectively, 
and is free of radial spokes. The Pad Chamber covers 
the fiducial pseudo-rapidity interval 2.05 < rj < 2.65 of the 
CERES spectrometer. It was added before the run period 
in 1995 as an external tracking device behind RICH-2 to 
assist the ring pattern recognition and reduce the fake-ring 
background in the high-multiplicity environment of Pb- 
Au collisions. By providing an absolute reference for the 
silicon and the RICH detectors, the Pad Chamber, with 
an angular resolution of about 0.6 mrad (in 9), proved a 
powerful tool in the geometrical inter-calibration of the 
detectors which helped to improve momentum resolution. 

The Pad Chamber is operated with a 90/10% Ar/C02 
mixture. The multi-wire anode is at equal distance (5 mm) 
to the upstream mesh cathode and the downstream pad 
cathode. Only about half the ionisation charges are col- 
lected during the 2/xs charge integration time of the pad 
readout electronics. The electronic avalanche produced by 
a charged particle traversing the Pad Chamber induces a 
signal in some of the « 29000 pads of the pad cathode. The 
pad size is the same as in RICH-2, and the pad readout 
electronics was adopted from RICH-2. 

The RICHes measure the change in local angle, while de- 
flection in the Pad Chamber is derived from the displacement 
relative to the distance from the vertex. 



2.6 The trigger 

A system of beam counters (BC) has been specifically de- 
veloped to meet the requirements of minimal mass expo- 
sure in the beam and target region and sufficient radiation 
hardness The trigger system is based on three small 
Cherenkov beam counters operating in air, one (BCl) 
about 60 m upstream, monitoring incident beam parti- 
cles, and one (BC2) about 6 m downstream of the tar- 
get detecting ions passing through the spectrometer. The 
third Cherenkov counter (BC3) registers each intact Pb 
ion downstream of the last target in order to veto the 
interaction trigger. To this purpose, Cherenkov photons 
(2700 per cm air) emitted in a narrow forward cone are 
refiected away from the axis onto a UV light guide by 
means of a tiny (6 mm diameter) aluminised mylar mir- 
ror on axis about 6 cm downstream of the target. The 
Cherenkov light is fed into a photo-multiplier just outside 
the spectrometer fFig l2.1|) . 

A plastic scintillator in front of the spectrometer (VC) 
is used to veto upstream interactions. The interaction trig- 
ger is defined as the logical AND of BCl and the veto of 
BC3 and VC, INT= BCl-VC -BCS. Centrality is selected 
with the multiplicity detector (MD), an array of 24 plastic 
scintillator paddles downstream of the RICH detectors at 
77 =2.9-4.7, the light output of which serves as a measure 
of the number of ionising particles that have passed. The 
centrality trigger is defined as INT-MD, and a hardware 
threshold is set at 100 mips. The accuracy of the trig- 
ger threshold and its stability over time is limited, mostly 
due to gain variations in the photo- multiplier tubes. In 
the off-line analysis, a precise multiplicity measurement is 
provided by the two SiDC's. The trigger selection corre- 
sponds roughly to the top 30% of the geometrical cross 
section fT^. A more precise calibration will be presented 
in sect. 13.81 



2.7 Data acquisition 

The on-line response of the detectors is shown in Fig. 12.21 
for a semi-central 158 GeV/n Pb-Au event recorded during 
data taking in 1996. 

The large amount of information per event contained 
in five highly granular detectors is collected by a fast data 
acquisition system. Data reduction is performed already at 
the hardware level when individual detectors are readout 
with zero suppression and pedestal subtraction. The data 
from the detectors are split into several separate readout 
chains which are processed in parallel for higher readout 
speed and transferred into special Memory modules where 
the data is compressed using Huffmann coding. 

A set of two (in 1995, three in 1996) CPU modules 
in VME technology is used to collect the data from the 
Memory modules for each event during the burst period 
(4.8 s) in a round-robin mode and writing it to tape during 
the intervals between bursts (14 s). As the write speed of 
the Digital Audio Tape (DAT) drives used is rather low, 
each CPU has three drives connected to it, allowing for 
an effective aggregated write speed of about 1.5 MByte/s 
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Fig. 2.2. Event displays for RICH-1 (top left) and RICH-2 (below), SiDC-1 (top right) and Pad Chamber (below). 



per CPU (in 1995, about 8 MB/s in 1996). This allowed to 
record on average 550 events/burst of 40-45 kByte each in 
1995; the upgraded DAQ together with optimised detec- 
tor settings and therefore smaller event sizes (30-35 kB) 
recorded about 1000 events/burst in 1996. 



2.8 Data taking in 1995 and 1996 

CERES/NA45 had data taking runs of 9 days in fall 1995 
and of 27 days in faU 1996 at the CERN SPS with a 
158 GeV/n Pb beam on Au targets. The average beam 
intensity in both years was about Ix 10^ ions per burst 
of duration 4.8 s. 

Because of the small target and beam dimensions, a 
readjustment of the beam position was one of the regular 
shift duties. To keep the spectrometer efficiency at the de- 
signed level, the gains in the UV detectors were monitored 
continuously and held within limits of about 30% at 2-10^ 
by adjusting high voltage upon changes of atmospheric 
pressure. 



Collision events were selected with the interaction trig- 
ger threshold set on 100 mips equivalent in the multiplicity 
array. The trigger contained an admixture of downstream 
interactions on the level of 15% of the target interactions. 
These were discarded in the off-line analysis on account 
of the silicon-drift track multiplicity. The latter served for 
more accurate centrality definition and also revealed some 
difference in the effective calibration of the multiplicity ar- 
ray between the two runs. 



In 1995 8.5 million events were collected with aver- 
age multiphcity {dNch/dr])= 220 corresponding to the top 
33% of the geometrical cross section. In 1996, the total 
sample was 42 million events of average charged particle 
multiplicity 250, or 26% of the top geometrical cross sec- 
tion. The multiplicity refers to the number of tracks in the 
SiDC's and is averaged over the pseudo-rapidity range 77= 
2-3. In the middle of the runs the polarity of the magnetic 
field was switched. 
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2.9 Instrumental means of coping with background 

We shortly review here the instrumental means by which 
CERES recovers a weak signal of low-mass electron pairs 
from high levels of combinatorial background. 

Approximate 'hadron blindness' is achieved by using 
two Ring Imaging Chcrcnkov (RICH) detectors with a 
high threshold ^th ~ 32. While electrons with momenta 
above 16 MeV/c produce Cherenkov light, pions overcome 
the threshold only at 4.5 GeV/c. More than 95% of aU 
charged hadrons pass without producing Cherenkov light. 

The radiation length within the spectrometer accep- 
tance has been kept at X/Xq 1 %. This is the level where 
the number of conversions is about equal to the number 
of Dalitz pairs. It is the result of persistent efforts in the 
design of all spectrometer components to reduce the de- 
tector materials, among which the thin mirror of RICH-1 
and the segmented target are the most important. 

Our physics sample are electron pairs with mass above 
200 MeV/c^. Below 200 MeV/c^, photon conversions and 
tt" Dalitz decays shoot up in yield which diminishes the 
sensitivity to interesting physics. In pursuit of the goal 
to recognise soft pairs of conversions and TTQ-Dalitz pairs 
with highest possible efficiency, the CERES spectrometer 
provides two powerful handles: the SiDC doublet detects 
close tracks by double-dE/dx response in pulse height, 
or by resolved close hits. The fact that the radiator of 
RICH-1 is free of magnetic field allows to see electron pairs 
of small opening angle undeflected, i.e. as close, resolved 
Cherenkov rings, or rings with a larger number of photons 
when the two tracks are separated by less than ~8 mrad. 
Soft electron tracks from conversions and 7r°-Dalitz pairs 
are strongly deflected by the magnetic field between the 
two RICHes, electrons and positrons in opposite sense. By 
setting an upper limit to the azimuthal deflection between 
RICH-1 and RICH-2, or between SiDC and Pad Chamber, 
a cut on track pt is implemented which is one of the most 
effective measures to reject background. 

To achieve a high efficiency for reconstruction of soft 
pairs, all detectors have full azimuthal coverage and the 
vetoing detectors have a slightly larger (minimum 1.9 < 
77 < 2.8) fiducial rapidity acceptance than the detectors 
after the magnetic field (2.15 < rj < 2.6). The more subtle 
details of the rejection strategy are discussed in sect. 13 



3 Data analysis 
3.1 Overview 



In this overview we sketch the strategy of the data analy- 
sis up to the pairing level when identified electron tracks 
in a given event are selected for combination to pairs. The 
presentation is mainly based on data and only rarely refers 
to Monte-Carlo simulations. The latter are, however, im- 
plicit in the choice of various quality and rejection cuts 
and will be treated in sect.^ 




100 200 300 400 500 
0ee(mrad) 



Fig. 3.1. Opening-angle and single-track pt distributions of 
7r°-Dalitz pairs (full histograms) and ui mesons (dashed his- 
tograms) in the CERES acceptance; Monte-Carlo simulations. 
Standard cuts (vertical lines) are not applied. 



3.1.1 The analysis strategy 

The invariant mass squared of the pair is given by the 
squared sum of the electron 4-momenta, 



Wee = (Pe+ + Pe 



y =2p,+ pe- (l-COs6)ee). (3.1) 



For the standard single-electron cut pt > 200 MeV/c used 
in the data analysis, the dynamic range of the electron mo- 
menta is somewhat restricted so that the dynamic range 
of the invariant mass is largely determined by that of the 
laboratory opening angle 0ee between the electron tracks. 

Signal electron pairs with m > 200 MeV/c^ have open- 
ing angles considerably larger than the asymptotic Cheren- 
kov ring radius of 30 mrad; the massive pairs from p, u, 
and 4> decays are opened about ten times wider (Fig. 13.11 
upper panel) . Our operational definition of signal or open 
pairs includes an opening angle cut 0ee > 35 mrad.^^ For 
77— Dalitz decays, the mean opening angle is about 60 mrad. 
Conversion and 7r°-Dalitz pairs have small masses^^ and 

It reduces the number of Dalitz pairs. As an important 
side effect, the cut enforces more uniform track distributions 
in polar angle of pairs from different sources. 

Photon conversion pairs are nearly massless but may ac- 
quire a small apparent mass by multiple scattering. 
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Fig. 3.2. Left: A 7r°-Dalitz decay, 7r° — > e"'"e~7, and a V track 
indicating a conversion pair, 7 — > e'^e~, show up in RICH-1 
as close doublets or single rings, followed in RICH-2 by two 
separate rings where one is missing (dashed), due to detector 
inefficiency or deflection out of the acceptance. Right: The re- 
maining single tracks are combined to a fake open pair which 
looks like a signal pair, opened already in RICH-1. 

average opening angles below 2 mrad and 20 mrad (rms), 
respectively. The sample of pairs with masses below 
200 MeV/c^ and 0ee > 35 mrad will be referred to as the 
Dalitz sample and is also used for checks on reconstruction 
efficiency and absolute yields. 

The pt distributions of electron tracks from conversions 
and 7r°-Dalitz pairs are steeper than those of open signal 
pairs. This feature provides the only rejection handle at 
the track level, albeit a very powerful one: the 200 MeV/c 
cut on track pt reduces close-pair tracks much stronger 
than signal tracks as seen in Fig. 13.11 lower panel. 

The most severe problem of the experiment is the enor- 
mous combinatorial background. We do not know which 
electrons belong to a pair and therefore we accept combi- 
nations of all tracks that qualify. However, when we find 
a pair with m < 200 MeV/c, its tracks are excluded from 
further pairing. Because the S/B ratio for these pairs is 
usually very good, we can declare a fully reconstructed 
Dalitz or conversion pair with good confidence. 

The background arises whenever low-mass pairs are 
only partially reconstructed and the remaining tracks are 
combined, as visualised in Fig. 13.21 Clearly, combinatorial 
pairs cannot be distinguished from genuine signal pairs 
and contribute to the entire mass range of interest, exceed- 
ing the signal by three orders of magnitude if all low-mass 
pairs would contribute. The final pair sample at masses 
above 200 MeV/c^ is still only about 10% of the resid- 
ual level of combinatorial background pairs. This is why 
already small inefficiencies in reconstruction of soft pairs, 
acceptance losses, etc., give rise to large relative levels of 
combinatorial background. Multiple scattering of conver- 
sions and 7r°-Dalitz tracks is a further important source 
of losses. 

A mild pt > 50 MeV/c cut is applied during the pro- 
duction stage (see below) to limit the search area in az- 



imuth. Many of the conversion and Dalitz pairs with only 
one leg above 200 MeV/c can be reconstructed and later 
rejected this way. Such pairs are ten times more numerous 
than pairs with both electrons above 200 MeV/c, so that 
many stiff electron tracks are taken out before entering the 
pairing stage. Applying the strong pt cut before the filter 
would have kept those tracks in the combinatorics. Tracks 
attributed to pairs with opening angles Oee < 35 mrad - 
the logical complement of the opening-angle cut for signal 
pairs - are marked to be excluded from further pairing. 

The deflection of electrons and positrons by the mag- 
netic field between the RICHes provides the unique search 
pattern of V-tracks, i.e. one ring in RlCH-1, possibly some- 
what blurred as it contains UV photons from two close 
tracks, and two separated rings in RICH-2 at about the 
same polar angle (see also Fig. 13. 2f) . 

Rejection of combinatorial background is optimised by 
tuning various cuts. Since high rejection power and high 
signal efficiency are competing requirements, an appropri- 
ate measure of signal quality is required. The optimisation 
has to be kept rigorously free of bias. A critical discussion 
of this important issue is given in sects. and 13.51 

3.1.2 The analysis stages 

Without a higher-level trigger in the CERES Pb-beam 
experiments, the actual data volume is much larger than 
that of interesting events, and a primary data reduction, 
the 'production' stage, is called for. Implementation of an 
effective production filter requires an accurate geometrical 
inter-calibration of all detectors. The availability of the 
Pad Chamber since 1995 allowed precise local tuning of 
the entire spectrometer using samples of high-pt pions j75j . 

During production, the full analysis chain for electron 
track reconstruction is at work, albeit under loose qual- 
ity criteria. Events which contain at least two electron 
tracks with pt > 50 MeV/c are stored in a database for 
further processing. As millions of events have to be pro- 
cessed with sophisticated pattern recognition algorithms, 
the production is time consuming. The 1995 production 
on the CS2 parallel computer at the CERN CN division 
with 32 SUN-SPARC2 processors took about 10 weeks; 
the 1996 data were preprocessed on a PC farm in several 
turns with readjusted production filters during 4 months. 
The final data analysis mainly deals with the optimisation 
of the pair sample. 

3.2 Reconstruction of electron tracks 

The reconstruction of electron tracks in the present analy- 
sis takes full advantage of the external tracking detectors, 
the doublet of SiDC's before the RICH spectrometer and 
the Pad Chamber after it. 



3.2.1 Coordinate systems 

Raw-data detector coordinates 
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In RICH detectors and in the Pad Chamber, hits are 
encoded as pad amplitudes in a two-dimensional mesh of 
{x,y) coordinates of the 'pad plane'. The natural unit 
is the pad size of 2.74 mm in RICH-1, and 7.62 mm in 
RICH-2 and the Pad Chamber, which is used up to the 
reconstruction of rings and ring centres in the RICHes, 
and of hits in the Pad Chamber. 

For the SiDC's, the symmetry of the radial drift field 
and the circular ring of anodes is maintained at the raw 
data level: the intrinsic coordinates are anode numbers 
(0-359) for the azimuth location and time bins (0-255) of 
20 ns for the drift time. 

Local detector coordinates 

RICH detectors measure angles of particle trajecto- 
ries that connect centres of Cherenkov rings with the ver- 
tex point, the units are A0 — As/ f, expressed by pad 
size As and focal length /; these are A6= 2.18 mrad and 
1.82 mrad in RICH-1 and RICH-2, respectively. The lo- 
cal coordinates in paraxial approximation are expressed 
by the tangents of the polar angle 6. Because of spherical 
aberration by the mirrors, the expression using 9 instead 
turns out to be a much better approximation, 

X = f 9 cos(f> + Xo, y — f 9 sincj) + Uo, (3.2) 

where Xq , ijo are the coordinates of the origin of the pad 
plane. In the off-hne analysis, the angles 9 and (j) are eval- 
uated by ray tracing using look-up tables to correct for 
local modulations in focal length. 

Hits in the silicon-drift detectors are given by their ra- 
dial position r (referring to the centre of the wafer) and 
azimuth angle (f> which is defined as in the global labora- 
tory system. The transformation from drift time to radius 
requires knowledge of the electron drift velocity with pos- 
sible spatial as well as temporal changes. 

Laboratory coordinates 

Once track segments of several detectors are to be 
joined, each detector is put into a three-dimensional global 
laboratory system with its own geometric calibration pa- 
rameters, i.e. small (x, y)-shifts of detector axes away from 
the optical axis, rotations around z, tilts, etc. This global 
laboratory system uses left-handed Cartesian coordinates 
with the z-axis along the beam and the x- and y-axes 
pointing to the right and upward, respectively, when look- 
ing with the beam. The origin is in the centre of SiDC- 
1. As we do not measure space points within the short 
magnetic field, it is convenient to work with straight tra- 
jectories all along from vertex to the Pad Chamber, and 
store the particle's momentum and charge sign derived 
from magnitude and sign of the azimuthal deflection. 

Event coordinates 

Eventually, in the global event coordinate system, the 
event vertex {xy, yvi ^v) is taken as the origin, and tracks 
are described by polar coordinates 9, 4> at the vertex. 

3.2.2 Hit finding in the SiDC's 

The 360 readout channels of anodes and the 256 time 
slices, or time bins, sampled by the FADC, span a matrix 
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each cell of which is assigned a 6-bit non-linear raw-data 
amplitude. The amplitudes are linearised and pedestals 
are subtracted. 

The data field is then searched for contiguous regions of 
cells with nonzero amplitude. Each cell is attributed to one 
such cluster. All cells of fixed anode number within a clus- 
ter form a time sequence of amplitudes called a pulse. The 
signal of a minimum-ionising particle produces pulses of 
typically 5 time bins and is spread over 2.2 anodes on aver- 
age. Pulses of less than 3 cells above a hardware threshold 
are discarded. Examples of such pulses on a few neigh- 
bouring anodes are displayed in Fig. 13.31 Due to comple- 
mentary signal transmission, the pulses are free of pickup 
which plagued previous runs. 

The centres of gravity of the pulses are calculated by 
Gaussian regression, or by Gauss fits taking the known, 
drift time dependent widths from a table. After time t, the 
drifting electron cloud arriving at the anodes has devel- 
oped a time-spread due to diffusion of cr"' = yJID t/vdrift- 
Here Z? ~ 35x 10^'*mm^//xs is the electron diffusion con- 
stant in silicon and Vdri ft denotes the drift velocity which 
varied in the range 6.0 - 8.5 mm//xs depending on the volt- 
age setting. The charge pulse is folded with the quasi- 
Gaussian response of the shaper, ^Yia,t the time 

spreads add in quadrature. The shaping introduces a bal- 
listic deficit in amplitude 6 — g-^haper j ^in -^^j^jf.}^ jg cor- 
rected for. Pulse heights saturated in the peak cells are 
approximately reconstructed. 

The stop-pulse correction mentioned in sect. 12.21 re- 
moves the random phase jitter. The conversion from drift 
time to drift distance is presented in sect. 13.21 




50 100 150 200 250 
time bin (20ns) 



Fig. 3.3. A collection of pulses on adjacent anodes. 1996 data. 
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Pulse trains on neighbouring anodes of the same clus- 
ter are merged into a hit if centres of gravity differ by 
less than one time bin. The hit coordinate in azimuthal 
direction is calculated as the centre of gravity of the con- 
tributing pulses, weighted by their peak amplitudes. The 
hit amplitude is the sum of the pulse amplitudes. 

3.2.3 Hit finding in RICHes and Pad Chamber 

Signals of the RICH detectors and the Pad Chamber are 
read out from the checkerboard-like arrangement of pads 
which receive the charge amplified in multi-wire propor- 
tional chambers. The same hit-reconstruction algorithm is 
used for these detectors. 

Adjacent pads with amplitudes above a readout thresh- 
old are connected to clusters some of which are caused by 
background and electronic defects. Typical background 
clusters, like long and thin stripes from ionising parti- 
cles on oblique trajectories, or clusters with many pads 
in saturation, are removed with the help of various clean- 
ing algorithms. Remaining clusters are split into regions 
containing one local maximum and are identified as UV- 
photon hits in the RICHes or charged-particle hits in the 
Pad Chamber. Hit centres are calculated as the centres of 
gravity of the contributing pads. 

3.2.4 Ring candidates in the RICH detectors 

RICH detectors require an additional step of pattern recog- 
nition to search for rings with asymptotic radius produced 
by electrons _77, . A typical RICH-1 ring is shown in Fig. 13. 41 
The ring search is done using the Hough transformation [JH] 
from the data field of the pad plane into the 'Hough ar- 
ray' (or 'parameter space') which has the same dimension. 
A given cell is the 'image' of all data points (hits) on a 
circle around that cell, and its 'Hough amplitude' is the 
sum of the data amplitudes on that circle. A ring candi- 
date shows up as a peak in the 'Hough array' which is 
the higher the more photon hits, or illuminated pads, lie 
on the ring. In practice, a 'digital' Hough transformation 
is used: all pads with signal amplitudes above a defined 
threshold enter with unit weight into the sum amplitudes, 
irrespective of their amplitudes. 

The event display in Fig. 13.41 demonstrates that rings 
can be formed also by random arrangement of single pho- 
ton hits, some of them produced by pious near the Che- 
renkov threshold. The first step in ring pattern recogni- 
tion is a linear Hough transformation of the pad plane 
onto the parameter plane. Besides the real maximum in 
the Hough array, there are other local maxima connected 
to fake ring centres. As a counter measure, a second, non- 
linear Hough transformation is performed which assigns a 
relative weight to each cell in the Hough array such that 
each hit counts most for its most favourable ring-centre. 
This way fake rings are suppressed, as can be seen from 
Fig. 13.41 by a larger gap that has developed between the 
amplitudes of the real and the fake rings. The parame- 
ter to select real rings at this stage of the analysis is the 
amplitude after the second Hough transformation. 




Fig. 3.4. (a) A Cherenkov ring in the RICH-1 pad plane, 
indicating pad amplitudes, (b) First Hough transformation, 
(c) Second, non-linear Hough transformation. 

All surviving ring candidates are assigned their final 
centre coordinates by a robust estimation which is based 
on iterative re-weighted least-squares fits of circles with 
asymptotic radius to the hits. A second fit with variable 
ring radius should eliminate charged pions which have a 
non- asymptotic radius. The function minimised is a mod- 
ified where the fit-potential varies in a Gaussian way 
(instead of quadratically) with the distance to the mini- 
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Fig. 3.5. Drift time spectrum, reflecting the pseudo- 
rapidity density of charged particle production. Small struc- 
tures are from local field variations. One time bin — 20 ns, 
Vdrift ~ 9.4 /im/ns. SiDC-2. 

mum. For an extensive description of the fitting algorithm 
see Refs. |7SII7^ . Other parameters besides the amplitude 
in the Hough array, like the number of hits and the spread 
of the hits around a perfect circle, can be used for fake-ring 
suppression. 

Ring reconstruction efficiency is determined by Monte- 
Carlo simulations with a cut on the number of photon hits, 
assumed to be Poisson-distributed. At the stage described 
up to here, ring reconstruction efficiencies of about 85% 
are still confronted with many fake rings, about 10 to 20 
for one reconstructed ring that is real. It was in antici- 
pation of such alarming majority of fake rings that the 
collaboration decided in 1994 to implement full external 
tracking 

3.2.5 Calibration of the SiDC telescope 

To first approximation, the relation between the drift dis- 
tance and drift time is linear. From the known radial ex- 
tension AR of the active area and the total drift time At, 
the drift velocity can be calculated as 

Vdrift = AR/At, AR = Rmax — Rmin, (3.3) 

where Rmin is the inner edge of the active area and Rmax 
corresponds to the anode radius. The corresponding 
drift time is 

At — tmax imin^ (3-4) 

where tmin is the time corresponding to the shortest drift 
path (ionisation directly under the anode), and tmax the 
drift time of particles starting at Rmin- A typical hit dis- 
tribution as a function of drift time is displayed in Fig. 13. 51 

^"^ The anode radius is 32 mm and 42 mm for the 3-inch and 
4-inch detectors, respectively, Rmin is typically 10 mm and the 
total drift time 3 to 5 fis. 



By fitting the edges of the drift-time spectrum, the values 
of tmin and tmax ^rc determined. 

Following such preliminary calibration, the geomet- 
rical alignment of both detectors and the determination 
of the interaction vertex is performed. Other corrections, 
such as the stop-pulse correction, corrections on drift- 
velocity variations due to temperature changes, etc., are 
done while maintaining the reconstructed vertex position. 



3.2.6 Vertex reconstruction 



Once hits are reconstructed, the information of hit posi- 
tions in the laboratory coordinate system is used to com- 
bine hits to tracks and find the vertex position to which 
almost all tracks of a given event point to. The proce- 
dure to minimise the quadratic sum of hit mismatches 
between the two detectors and its iteration is extremely 
time consuming. We used therefore a robust vertex fitting 
approach 80 : all hits in SiDC-1 and SiDC-2 (typically 
more than 100) are combined to straight track segments 
and a weighted sum of their projected distances to the as- 
sumed vertex position is calculated. In the next iteration, 
this centre of gravity becomes a new starting value for the 
vertex position and each track segment gets a new weight 
according to its deviation from the mean value in the step 
before. After the position of the vertex is determined, its 
z-position is redefined as the exact position of the closest 
vertex disc. 

Figure I^Sl displays the density of reconstructed vertex 
positions along the beam. The peaks reveal the positions 
of the eight discs of the segmented target assembly with 
a spacing of 3.2 mm. On average, the reconstruction was 
done with 160 charged particles per event. The data shown 
in Fig. l3.6l was accumulated over 2.6 x 10^ Pb-Au events, or 
six hours, and demonstrates a certain long-term stability. 
The resolution^* of = 250 /zm is sufficient to identify 
the correct target disc without ambiguity^^ . 

More critical for tracking accuracy is the precision with 
which the vertex can be localised in the plane transverse 
to the beam, as there are no fixed points. To measure it 
we used stiff pion tracks with pt > 1.2 GeV/c to minimise 
multiple scattering. The scatter plot in Fig. 13.61 accumu- 
lates over many events the distance in x and y between the 
actual event vertex (determined from all tracks) and the 
point where a given pion track intersects the respective 
target disc. Projecting on the axes, we obtain the 1-dim 
distribution shown to the far right; a lateral resolution of 
28 /im (x and y) is derived. The transverse vertex res- 
olution is therefore ar = V^ctx — 40 fim. By choosing 
the vertex position in the x, y-plane as the origin of the 
event coordinate system, we account for event- by- event 
displacements within the diameter of the target. 



i.e. the standard deviation in the mean z position 
due to the unknown location of the vertex inside the 25 jj.m 
thick foil, there is an rms error of 25/vT2, about 7 /im. 
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Fig. 3.6. Reconstruction of interaction vertices, (a) Density of ^-positions with peaks corresponding to the eight target discs. 
(b,c) Vertex resolution in the x-y plane determined from high-pt pions. 1996 Pb-Au data. 



3.2.7 Charged particle tracks 



3.2.8 Straight pion tracks for calibration 



Silicon-drift track segments are constructed by connect- 
ing the vertex point to hits in SiDC-2 which lie within the 
fiducial acceptance. A track segment is accepted if there 
is at least one hit in SiDC-1 within a predefined window 
around the point of intersection; for more than one hit, 
the centre of gravity is taken. Once the interaction vertex 
and the SiDC track segments are reconstructed, the tra- 
jectories of charged particles are extrapolated downstream 
to the Pad Chamber. If a pad hit is found within a certain 
fiducial window, a track candidate is created. 

The sizes of the fiducial windows are expressed as mul- 
tiples of the rms widths of the corresponding matching 
distributions; usually these are 5 Gmatch during produc- 
tion and ?> (J match during the final analysis. The fiducial 
windows were taken momentum-dependent to reduce effi- 
ciency losses at low momentum due to multiple scattering. 

The matching of tracks in direction between detec- 
tors separated by the magnetic field requires special at- 
tention since such tracking window corresponds to a mo- 
mentum cut. To avoid loss of tracks by multiple scattering 
for large 4> deflection, the fiducial window in direction is 
opened to assume the shape of a butterfly. During the pro- 
duction stage, we use a tracking window in (/) of ± 0.6 rad. 
This corresponds to the pt-cut of 50 MeV/c mentioned al- 
ready. During off-line analysis, the size of the butterfly is 
approximated by a sector of 100 mrad in 0, corresponding 
to a lower momentum cut-off of 1 GeV/c, the standard pt 
cut of 200 MeV/c, and 3 mrad in direction. Matching 
distributions will be discussed in sect. 13.21 

The contribution of background hits was estimated by 
applying detector rotations. Since the hits in the SiDC's 
are highly correlated, the main background contribution 
comes from random combinations of SiDC track segments 
with hits in the Pad Chamber. By rotating these hits in a 
given event by a random angle with respect to the silicon 
detectors, the true physics signal is destroyed and only 
background tracks remain. 

By requiring that track elements from all five detectors 
match within three standard deviations of the detector 
resolutions combined with the rms spread due to multiple 
scattering, fake tracks are reduced to a negligible level. 



Straight tracks of high-momentum pions are an impor- 
tant tool for fine-tuning the spectrometer calibration l75] . 
Pion identification uses external tracking by the SiDC's 
and the Pad Chamber to predict ring centre positions in 
the RICHes. Pions have been selected under tight qual- 
ity criteria regarding matching, number of photons on 
Cherenkov rings, and clean environment around rings. 

To extend the SiDC track segments into the Pad Cham- 
ber, only a narrow window of ±30 mrad in azimuth is 
searched for a matching hit, variations in deflection be- 
ing very small. In polar direction, the Pad Chamber is 
searched merely over the matching window between the 
two detectors since multiple scattering is negligible. All 
hits found are candidates. From the coordinates before 
and after the field, the ring centres in the RICH detectors 
are predicted. The pointing to RICH-1, without deflec- 
tion, is unproblematic. With the momentum information 
derived from the azimuthal deflection between SiDC and 
Pad Chamber, the tracks are extrapolated into RICH-2, 
behind the field. 

Photon hits that fall into the vicinity of the predic- 
tors are collected and used as input for a robust fitting 
algorithm jZJj. To find centre and radius of the ring, the 
rms deviation of hit positions from a circle with radius 
R around the ring centre is minimised in a Gaussian fit 
potential with three free parameters. Once the radius R 
is found, the pion momentum is re-evaluated from the re- 
lation 

. m 

(3.5) 



1th; 



^l-iR/RoofY 



where i?oo stands for the asymptotic ring radius. 



3.2.9 Spectrometer calibration 

For the analysis of the 1995 and 1996 data every spectro- 
meter component was calibrated. Several successive steps 
were necessary to align all detectors and to derive appro- 
priate correction functions for the analysis chain. A first 
rough calibration was performed to adjust detector offsets 
and to correct small rotations and tilts. 
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Fig. 3.7. Vertex positions of the 8 target discs taken when 
cooling of the SiDC's was switched off (left), and after stabil- 
ising on the drift time spectrum (right). SiDC-2, 1995 data. 
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Fig. 3.8. Deviations in radial hit position with respect to 
reference of stiff n tracks; before calibration (left, full circles), 
after tuning the iid,i/t-parameter (left, open circles), and the 
tmin parameter (right). SiDC-1, 1995. 



Since electron mobility in silicon strongly depends on 
temperature^*^, the calibration is strongly affected by tem- 
perature variations. To keep the calibration stable over pe- 
riods of hours, we employed the simple and fast method by 
which the position of the upper edge of the drift time spec- 
trum is monitored to provide an 'online' drift-velocity sta- 
bilisation. The feedback procedure successfully stabilises 
the calibration as can be seen from the resulting stability 
of vertex positions in a set of test measurements taken 
while the cooling was switched off. Fig. 13.71 Under nor- 
mal running conditions, the temperature variations were 
typically below 0.5°C over 12 hours. 

On even longer time scales, drifts in the calibration of 
the SiDC's were prevented employing the fixed reference 
provided by the Pad Chamber using high pt high-statistics 
pion samples. The method registers the deviations of hits 
in either one of the SiDC's from the straight line that 
connects the event vertex with a selected hit in the Pad 
Chamber; it works fine in a pre-calibrated system. The 
effects of fine tuning the parameters Vdrift and tmin of 
SiDC-1 can be seen in Fig. EH 

Response from RICH-1 when given close scrutiny by 
the absolute reference grid of the Pad Chamber revealed 
slight deformations of the spherical mirrors which caused 
non-linearities in polar angle. Following the observation 
that the radius of curvature of mirror- 1 decreases towards 
the rim, application of a small ^-dependent variation in 
the focal length removed this problem. 



3.2.10 Matching quality 

Starting with an internal calibration of the SiDC vertex 
telescope, the remaining detectors are aligned to the cen- 
tre of SiDC-1. Internal consistency and quality of the read- 
justed calibration can be evaluated from the residual offset 
in centres of gravity, the widths of the matching distribu- 
tion close to the peak, and the amount of background 



The electron mobility jj. 
ture as /i oc T"^ *. 



Vdrift 



/E depends on tempera- 



in the tails, for all detector combinations. Examples of 
matching distributions are shown in Fig. 



At low momenta, the matching quality deteriorates 
because of multiple scattering in the material, and it is 
no longer determined by the detector resolutions alone. 
The widths of the matching distributions in 6, plotted for 
SiDC - Pad Chamber and RICH-1 - Pad Chamber shown 
in Fig. I3.1UI for measured and simulated spectra, exhibit 
very well the dominance of the multiple-scattering contri- 
bution increasing with 1 /p over the momentum-independent 
detector resolution towards small momenta. The simula- 
tions with the actual spectrometer parameters are shown 
by dashed lines. Final track selection is achieved by ap- 
propriate momentum-dependent cuts in the matching dis- 
tributions. 
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Fig. 3.9. The ^-matching distributions for various detector 
combinations determined with high-pt charged pions. 
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Fig. 3.10. Quality of track matches A9 of Pad Chamber to 
RICH-1 and to SiDC vs momentum. Shown are 1996 data and 
Monte-Carlo simulations. See text. 

3.3 Momentum analysis and mass resolution 

Since particles passing the SiDC's and RICH-1 have ex- 
perienced no magnetic field, these detectors are used as 
zero-deflcction reference. The deflections Acj) of charged 
particles of momentum p between the RICHes and be- 
tween the SiDC's and the Pad Chamber were given in 
sect. 12.41 The best resolution is obtained from all pos- 
sible detector combinations, weighted by the respective 
accuracy. The result of the simulation with the measured 
detector resolutions as input is shown in Fig. Kllll By de- 
tailed Monte Carlo simulation which realistically describes 
the individual detector resolutions and the quality of the 
matching distributions, the overall resolution is described 
by the function 



Ap/p = V(2.3%-p)2 + (3.5%)2, (3.6) 



Q. 

-a 



p(GeV) 

Fig. 3.11. Momentum resolution obtained with different de- 
tector combinations vs momentum; MC simulations. The solid 
line combines the information of all five detectors. 
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Fig. 3.12. Correlated SiDC amplitudes in 5 mrad summing 
window for sample enriched in photon conversions. Majority of 
particles have single tracks in both SiDC's (pions), or double 
amplitudes which signify electron pairs from photon-conversion 
in the target. 

with p in GeV/c. The momentum dependent part is due 
to detector resolutions, the constant part due to multiple 
scattering. The mass resolution at the p/uj is about 6%, 
at the (j> about 7%. 

3.4 Rejection of combinatorial background 

With a.pt cut of only 50 MeV/c at the production level, an 
important part of the rejection is already achieved: very 
soft tracks are removed, the tracks of all reconstructed 
conversion and Dalitz pairs are marked and taken out of 
the further analysis, including pairs with only one electron 
of pt > 200 MeV/c. Only then is the single-track pt-cut 
tightened from 50 to 200 MeV/c which reduces the com- 
binatorial pair background by about a factor of 10 while 
keeping 97 % of signal pairs. Tracking cuts are narrowed 
down from 5a to 3cr. Together with increased requirements 
in track and ring quality, this results in a drastic suppres- 
sion of fake tracks as mentioned before. 

The most powerful rejection tools on the pair level de- 
rive from the fact that target conversions and a large frac- 
tion of 7T° Dalitz decays produce close tracks. These are 
either registered as single hits of 'double-dE/dx' in the 
SiDC's when the separation is less than about 3 mrad, or 
as a 'double ring' in RICH-1 when the two rings coalesce 
into a single structure for separations below 8 mrad; or as 
more or less well resolved tracks, still considerably closer 
than the mean spacing of charged particle tracks in the 
SiDC's or of electron tracks in RICH-1. 

To maintain high rejection power also in the region be- 
tween perfect overlap and resolved double hits, we use a 
'summation' window in the SiDC's of 5 mrad in which hit 
amplitudes are summed up'^^. Figure r^.l2\ shows a two- 
dimensional plot of the hit response in SiDC-1 vs SiDC-2 
obtained this way. It demonstrates that with an appro- 
priate two-dimensional rejection-cut, photon conversions 

■^^ This allowed to relax demands on double-hit resolution and 
thereby to avoid excessive (and mysterious) 'hit splitting'. 
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in the target are efficiently rejected without cutting much 
into the Landau tail of single-hit distributions. RICH-1 
identifies close tracks of conversions and Dahtz pairs by 
the number of photon hits, which suffers severely from 
pile-up losses, however. Alternatively, we use directly the 
analog sum of the gain-corrected pad amplitudes to iden- 
tify 'double rings'. Photon conversions occurring in SiDC- 
1 differ in response from target conversions only by pro- 
ducing on average a single-hit signal in SiDC and can be 
recognised by such signature. 

To be specific, we use the following four rejection steps 
which evolve in complexity with the information gained 
along the trajectories: 

1. Tracks are rejected if the amplitudes of the 5 mrad 
summation window in both SiDC's are larger than the 
typical single hit response encoded by parameters 
Slfiigh and S2high- 

2. Conversions in SiDC-1 are rejected by the requirement 
that the amplitude in SiDC-1 is below Sliow and the 
amplitude in SiDC-2 above S2high, summed over a 7.5- 
mrad window; and by observation of a double ring in 
RICH-1. 

3. A track with a double ring in RICH-1 is rejected al- 
ready if Shigh is surpassed in only one of the SiDC's. 
Besides improving rejection efficiency in general, this 
cut rejects conversions in SiDC-2. 

4. Tracks are rejected if another ring in RICH-1 is found 
within a wider window of 35 mrad. To avoid excessive 
vetoing by accidental structures, the second ring is re- 
quired to be of high quality and to connect to a track 
segment in both SiDC's. 

A few comments are in order. The search in the SiDCs 
for hits near tracks cannot be extended to larger distances 
because of the increasing chance to find a pion, vetoing 
the signal. For the purpose of 'Dalitz rejection', the veto 
is better based on RICH-1 where the inspection area can 
be opened up to 35 mrad without much signal loss because 
the density of electron tracks is so much lower. The strat- 
egy described above is quite successful in rejecting conver- 
sions but even with the close-ring cut it is of limited power 
in rejecting 7r''-Dalitz pairs because of the larger opening 
angles involved. 

In the last step we also reject surviving charged pion 
tracks. Figure 15'. 131 is a scatter plot of track deflection Z\0 
as abscissa, and the sum of the reduced ring radii in both 
RICHes on the ordinate; this quantity only depends on 
the Lorentz 7 factor. Because of different mass compared 
to electrons, pions of same momentum (deflection) have a 
different 7 (radius). It is seen that pion tracks even with 
asymptotic radius are clearly distinguished from electron 
tracks by virtue of the much smaller deflection in the field, 
and they can be rejected by a two-dimensional cut in both 
RICH detectors on ring radius vs deflection. This cut was 
used only in the '96 analysis; it improves the signal-to- 
background ratio at large masses and has little impact at 
lower masses. 
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Fig. 3.13. Scatter plot of the sum of the reduced ring radii 
R/Ra^ in RICH-1 and RICH-2 over deflection in the field for 
reconstructed electron tracks with asymptotic radius Raa- Ver- 
tical stripes are from charged pion tracks of both charge signs. 
The pion rejection cut used is indicated. 



3.5 Optimising signal quality 

In tuning the parameters of quality and rejection cuts, 
we strictly avoided to optimise on the data itself. Such 
practice is known to increase the risk of statistical fluc- 
tuations producing spurious, misleading results. Rather, 
the signal efficiency was determined by overlaying Monte- 
Carlo-generated tracks chosen from the hadronic cocktail 
on real events and measuring the rate of successful recon- 
structions using the quality and rejection cuts as in the 
current step of the data analysis. The signal efficiency may 
also be monitored and optimised using a sample of fully re- 
constructed Dalitz pairs The rejection steps were 
also tuned by Monte-Carlo techniques with generated con- 
versions and 7r°-Dalitz decays. Such a Monte-Carlo sim- 
ulation of the background is quite demanding and was 
only achieved in the 1996 data analysis 56 . In the 1995 
analysis, the background was taken from the data sample 
itself. 

By these procedures, each analysis cut is arbitrated 
objectively according to how much in background rejec- 
tion is gained for how little loss in signal efficiency. A clear 
understanding of what is meant by an 'optimal' balance 
is provided by the 'equivalent background-free signal'^^ 



For a given rejection cut which reduces the signal from S 
to S' = £ S and rejects a fraction {B — B')/B = r of the 



It is equal to the number of signal pairs of the same statis- 
tical significance if there would be no background. 
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background, the ratio 



S" 

^ eq 
Sea 



1 



(3.8) 



should be maximised. 

In summary, the rejection steps together reduce the 
number of background pairs by factors of 15 and 12 for 
mass above and below 200 MeV/c^, respectively, while 
signal pairs are reduced by only 20%. 



3.6 Pairing and subtraction of combinatorial 
background 

The combinatorial background B of unlike-sign pairs can 
be accounted for by the like-sign pair sample because by 
lack of any physics source the latter is purely combina- 
toric, 

(3.9) 



B = 2 yjN++N- 



Here, N++,N denote the numbers of e+e+ and e~e~ 

pairs, respectively. The relation is exact when charge sym- 
metry is fulfilled. The number of signal pairs S is obtained 
by subtracting the number of background pairs from the 
number of unlike-sign pairs, 



S = - B. 
The variance in S is approximately^"^ 



N+- 



B. 



(3.10) 



(3.11) 



Typically, the open pair signal S accounts for only a small 
fraction of our total unlike-sign pair sample for masses 
above 200 MeV/c^. 

The asymmetry induced by Compton electrons and 
decays is below 1% and would alter the factor of 
2 in eqn. (3.12) by less than 1 x 10^"^. It is worthwhile 
to note that eqn. (3.12) remains a very good assumption 
also in case of a small imbalance in the number of tracks 
of positive and negative charge, be it due to different ac- 
ceptance or reconstruction efficiencies, or to a generic im- 
balance j55) . For an asymmetry of 10% in the number of 
positive and negative tracks, the factor 2 in the above ex- 
pression would decrease by 2% only. 

The asserted symmetry has been tested by generating 
purely combinatorial pairs of either type, like-sign and 
unlike-sign, from data by 'reversing' the rejection cuts. 
The recipe is the following: prepare a clean sample of con- 
versions and Dalitz decays and recombine all tracks such 
that none meets its original partner. This sample is void 
of signal pairs and in every respect resembles the experi- 
mental pair background. The mass spectra of both com- 
ponents of this artificial background sample are compared 
with each other in Fig. 13.141 No significant charge asym- 
metry is visible between the mass spectra of the like-sign 



Combinatorial pairs are non-Poissonian, their variance ex- 
ceeds the mean. Because the density is only of order 10"'' per 
event, the deviation from eqn. (3.14), however, is below 5%. 
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Fig. 3.14. Comparison of mass spectra of strictly combina- 
torial pairs generated from clean data samples of conversions 
and Dalitzes (see text). Neither in form nor in absolute yield 
does a significant charge asymmetry become visible. 



and the unlike-sign pair background, neither in shape nor 
in absolute yield. Over the full mass range, the ratio of 
the integral yields is 



>;+e-/(i^e+e+ +>;-e-) = 0.996 ± 0.005, 



(3.12) 



which excludes an asymmetry larger than 1% at 90% con- 
fidence. We shall return to this topic in sect. 16. 5( 

Let us take an alternative route of estimating combi- 
natorial background. A philosophically correct way is to 
generate unlike-sign pairs by combining strictly uncorre- 
lated tracks chosen from different events under the same 
kinematic constraints as applied to the data. This is com- 
monly referred to as the mixed event method. Practically, 
efficiency losses that occur in the true data for close tracks 
have not been accounted for. A comparison of combinato- 
rial mass spectra obtained by the two methods from the 
full 1996 data set is displayed in Fig. 13.151 The coarse 
comparison on a log scale (top panel) does not show de- 
viations. On a finer linear scale (bottom panel), the ratio 
indicates that the like-sign pairs loose up to 30% compared 
to the mixed-event generated unlike-sign pairs - which is, 
however, due to the mentioned loss in detection efficiency 
at close track separation; an effect not present in the cur- 
rent realization of the mixed-event method. Yet, the ana- 
lysis allows to conclude that above m « 300MeV/c^ the 
ratio stays at unity within a band of less than ±1.5%. 

Measured mass spectra of like-sign and unlike-sign pairs 
are shown in Fig. 13.161 For comparison, a like-sign pair 
background has also been generated by random pairing 
of electron tracks with momenta, polar angles and trans- 
verse momenta chosen at random from measured distribu- 
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Fig. 3.15. Comparison of combinatorial background spectra 
from same-event like-sign pairs and, alternatively, from unlike- 
sign pairs of different events. Top: Comparison on log scale. 
Bottom: ratio 'same'/'mixed' on linear scale. 1996 data. 



tions. The result of such simulation is shown in the figure 
and compares very well with the measured like-sign mass 
spectrum. The simulated opening-angle distribution (not 
shown) also agrees well with the measured distribution. 

Figs. 13.141 13.151 and 13.161 display the typical shape of 
the combinatorial mass spectrum for the applied cut on 
track pt ■ This kinematical condition produces a falloff to- 
wards low masses and a broad peak at about twice the 
minimum track pt . The peculiar shape of the combinato- 
rial mass spectrum peaking around mass of 500 MeV/c^ — 
where the enhancement dominates — raises doubts whether 
the combinatorial background has been correctly assessed 
and subtracted. We return to this issue in sect. 16.51 to as- 
sert an unbiased background handling. 
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Fig. 3.16. Unlike-sign (open squares) and like-sign mass spec- 
tra (full circles) after all rejection cuts. Like-sign combinato- 
rial background shown as smooth curve is obtained by random 
pairing of tracks with the same kinematical composition as in 
data. The drop in yield below 400 MeV/c^ originates from the 
single track pt-cut of 200 MeV/c. 1996 data. 



3.7 Spectra and statistical errors 

The raw signal counts and their relative statistical errors 
are obtained by subtracting the combinatorial background 
channel by channel from the spectrum of unlike-sign pairs, 



Si = ini,+- - bi), ai/si = ^n.i,+_ + bi/{ni^+^ - hi). 

(3.13) 

The counts per channel add up to the total counts known 
already from eqns. (3.12) and (3.13), 



B. 



(3.14) 



Pair yields per event, corrected for reconstruction effi- 
ciency and normalised to {Nch), are written symbolically 



(3.15) 



with g a scaling factor and 1/e for efficiency correction, 
possibly channel dependent. Relative statistical errors are 
determined by the raw data counts. 



(3.16) 



To reduce staggering in signal counts due to channel 
fluctuations in the background spectrum, 'smoothed' com- 
binatorial pair spectra free of bin-to-bin fluctuations have 
been subtracted. The method allows, moreover, to limit 
the fluctuations in data points of background-subtracted 
spectra to truly statistically independent errors; it is dis- 
cussed in the following. 

The Monte-Carlo generation of background spectra is 
described in the previous section. The subtraction from 
the unlike-sign pair spectrum is done bin by bin. 



i , + - 



-hB. 



(3.17) 



The fractions fi add up to 1. The bin-to-bin errors 



sjnooth 



(3.18) 



in our case of B / S 3> 1 amount to a fraction quite accu- 
rately of l/\/2 of the full error of (3.19). The finite sample 
error in B, to which the smooth background spectrum is 
normalised, was neglected so far. It contributes to each 
channel a share of the same magnitude as the bin-to-bin 
error itself^^. 

Since the 'normalisation' errors common to all chan- 
nels differ in size, according to the magnitude of fi in 
eqn.(3.20), they cannot be represented as an error in scale. 
However, there is no need to show them separately; we 
only have to keep in mind that the data points have nor- 
malisation errors which are of the same size as the sta- 
tistical bin-to-bin errors which will be displayed, but are 
tightly correlated within the entire spectrum. 



because these errors add in quadrature. 
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Fig. 3.17. Distribution of centrality triggers for the two data 
samples of 1995 and 1996. Centrality is expressed by the num- 
ber of reconstructed charged particle tracks in the two SiDC's 
in the range 2 < y < 3, corrected for pile-up and efficiency. 
The selection corresponds to the most central 33% and 26% of 
the total inelastic cross section for 1995 and 1996, respectively. 

3.8 Centrality determination 



o data 

— UrQMD minbias 

- - UrQMD mix to 26% 




(2<il<3) 

Fig. 3.18. UrQMD: minimum-bias N^h distribution scaled 
to fit the upper corner (dashed histogram), and distribution 
of three centralities composed to describe the trigger profile at 
the lower corner (solid histogram); compared to the triggered 
Nch distribution (string of small circles), see text. 1996 data. 



The measurement of absolute cross sections is hampered 
by the fact that the number of particles that actually pass 
the segmented target is not precisely defined; this is be- 
cause the diameter of the Au disks is not much larger than 
the beam diameter. We decided therefore to calibrate the 
centrality of collisions by the shape of the pseudo-rapidity 
density of charged particles as measured by the SiDC's. 
Once the measured Nch distribution is corrected for var- 
ious instrumental distortions, it is used to calibrate the 
charge density obtained from minimum-bias UrQMD cal- 
culations IHJ. These still differ from the measured distri- 
bution at the \ow-Nch tail due to the non-ideal response of 
the interaction trigger. The final step consists in describ- 
ing the triggered Nch distribution by a linear combination 
of UrQMD calculations belonging to different impact pa- 
rameters. These steps are described below. 

3.8.1 Charged-particle density 

Charged-particle density in the range 2 < 77 < 3 is derived 
off-line from the number of tracks that emerge from the 
event vertex and intersect the two SiDC's sufficiently close 
to reconstructed hits. The determination of the track re- 
construction efficiency is a rather complex problem. To 
the percent accuracy required, it is influenced not only by 
hardware imperfections (i.e. dead anodes, electronic noise, 
pulse shape) , but also by the quality of hit and track re- 
construction. In particular, 'pile-up' effects due to finite 
double-hit resolution and artificial hit splitting are the 
most important effects to be corrected for. To that goal, 
we simulate charged particle tracks in the SiDC's gen- 
erated from UrQMD events. The analysis of the Monte- 
Carlo events uses the standard data analysis software, be- 
sides larger matching windows of 5 (Jmatch- Details are 
described in sect. ^ 



The functional dependence of the inverse of the recon- 
struction efficiency e on the measured number of charged 
particles is of the type 

{Nch,measured^ — ^ch,true / -^ch.measured — 
a ~\~ b Nch,7neasured- (3.19) 

The intercept at Nch.measured = corresponds to the in- 
verse 'static' efficiency and includes not only losses 
due to dead anodes, but also gains from artificial hit split- 
ting. It amounted to 97% and 93% in 1995 and 1996, re- 
spectively. The parameter b is about (6-7)xl0~^ which 
corresponds to about 10% relative losses, considerably more 
than expected from pile-up alone. 

The corrected multiplicity distribution of the trigger- 
selection is shown in Fig. 13.171 In the middle part and at 
the high-A^ch side, the shape of the distribution is sup- 
posed to be an undistorted image of partial cross sections, 
but not so at low Nch where the trigger profile of the 
MD with threshold set at 100 mips becomes visible (see 
sect. I2.6|l . The thresholds of the centrality trigger corre- 
sponded roughly to 30% for '96 and 35% of the total in- 
elastic cross section for '95. 



3.8.2 Calibrating the centrality using UrQMD calculations 

A more precise calibration of the centralities is described 
in the following. Minimum-bias distributions obtained by 
UrQMD (version 1.3) have been scaled slightly to make 
their upper corners coincide with the measured trigger 
distributions. From the UrQMD minimum-bias collisions 
several centrality classes are sorted out according to cer- 
tain impact parameter ranges. To reproduce the shape of 
the trigger profile at the lower edge of the Nch distribution. 
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a suitable linear combination of the corresponding Nch 
distributions is used. The resulting distributions are dis- 
played in Fig. 13.181 The total inelastic cross section of 
6.94 barn, and the trigger fractions were calculated with a 
geometrical overlap model. The resulting x = a juinei are: 

a;('95) = 0.6 x 35%-}-0.2 x 33%-l-0.2 x 28% = 33% (3.20) 

x('96) = 0.16 X 35% + 0.15 x 30% + 0.69 x 23% = 26%. 

(3.21) 

The systematic error from variation in the linear combina- 
tions alone is estimated as ±1.5% absolute. For the unified 
data, an average centrality of 28% of the top geometrical 
cross section has been adopted with an estimated uncer- 
tainty of about ±2%. 



4 Monte Carlo simulation 

4.1 Detector simulation 

The spectrometer with all material in proper geometry has 
been implemented in the GEANT 83 detector simulation 
package version 3.15. The present version of GEANT was 
not able to describe the number of photons per electron 
ring correctly, which forced us to write our own function. 

The detector response is simulated by taking into ac- 
count the particular signal generation including all known 
effects that influence position, width and amplitude. In 
the SiDCs these include the charge division among adja- 
cent anodes, local variations in drift-velocity over radius, 
diffusion along the drift path in radial and orthogonal di- 
rection, anode-wise gain variations including dead anodes, 
electronic noise from the readout, and digitisation errors. 
In the RICH detectors we account for chromatic aberra- 
tion and mirror quality. In the UV detectors and the Pad 
Chamber, simulations include the transverse diffusion in 
the conversion zone and the first amplifying gap, the digi- 
tisation effects in the last multi-wire amplification, and 
the noise of the readout electronics. 

Most crucial but difficult to achieve is a realistic sim- 
ulation of the background in which the simulated signal 
is to be embedded. The background is caused by charged 
particles, gamma rays and electro-magnetic radiation pro- 
duced directly or indirectly by the collision, ^-electrons, 
and particles passing the UV-detectors. In addition, the 
entire spectrometer or parts thereof might deviate from 
response linear with produced signal charges due to tran- 
sient saturation effects of various kinds. No sufficiently 
reliable simulation of the background was possible up to 
now. We use real data events into which the simulated 
Monte Carlo (MC) signals are embedded. This is done 
by adding the MC amplitude for each single cell, or pad, 
on top of the data amplitude. The systematic error in- 
troduced by increasing the detector occupancy has been 
estimated to be negligible. 

The MC signal was reconstructed using algorithms iden- 
tical with those of the data analysis. By comparing the sig- 
nals of reconstructed MC tracks and reconstructed data 
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Fig. 4.1. Comparison of reconstructed Cherenkov rings in 
Monte-Carlo simulations (stars) and in data analysis (his- 
togram) regarding the number of hits per ring (top row), the 
sum amplitude (middle row) , and the amplitude of the second 
Hough transformation (bottom row). 



tracks we fine-tuned the parameters of the detector sim- 
ulations to achieve the best description of the data. The 
characteristics employed for the comparison were the am- 
plitude distributions of hits, the widths, the numbers of 
pads/time-bins/anodes contributing to a hit, the num- 
ber of hits belonging to a ring and the local variations 
of all these parameters over the detectors. An example of 
the agreement achieved between MC simulation and data 
analysis in reconstruction of Cherenkov rings is shown in 
Fig. EI] 

It is crucial to model quantitatively the response to 
compound signal patterns, like accuracy in reconstruction 
of hit positions or ring centre positions, or the number dis- 
tributions of reconstructed photons per ring. Quantities of 
that kind depend also on characteristics of the physics in- 
put like opening-angle distributions of pairs in the sample 
and the polar distribution of electron tracks. They are 
therefore best checked on the data itself. For that pur- 
pose a special sample enriched in Tr'^-Dalitz decays was 
extracted from the data by searching for pairs with open- 
ing angles less than 50 mrad 0^1 • Such pair sample has a 
typical signal-to-background ratio of order one. 

The properties of this sample compare well with the 
corresponding properties of simulated Tr'^-Dalitz decays in 
ring-parameters and efficiency. There is but one exception: 
the measured ring centre resolution of the first RICH de- 
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tector is worse by a factor of 1.4 compared to the MC 
simulation, and also deviates from the expected resolution 
for single photon hits. For lack of understanding, wc in- 
troduce an ad hoc Gaussian random smearing of 0.9 mrad 
into the simulation to meet the measured value. The influ- 
ence on final momentum resolution is very weak, since the 
azimuthal track position measured in the SiDC's is any- 
how more accurate and outweighs the RICH contribution 
to tracking. 

4.2 Determination of reconstruction efficiency 

To compare the number of reconstructed electron pairs 
to any physics-based expectation, an absolute normalisa- 
tion is required. To correct spectra for reconstruction ef- 
ficiency^^, measured yields are multiplied by the inverse 
pair reconstruction efficiency averaged over the full sample 
of pair candidates. We may call this a correction 'on a sta- 
tistical basis', in contrast to an 'event-by-event correction' 
where the invariant-mass or pf^ spectra are incremented 
with the inverse efficiency of that event, or of particular 
track candidates. 

One way to determine the reconstruction efficiency is 
to measure the number of pairs from a well-defined physics 
sample, and adjusting it to the expected number. We first 
explain this 'Dalitz method' and then turn to the alterna- 
tive Overlay Monte-Carlo technique. 



Note that only the cocktail part factorises into averages 
of electron pairs and charged particles^^ while the exper- 
imental average does not, since Nch/N^f, depends, via re- 
construction efficiency on Nch- 

The recipe that follows is to average the ratio of the 
number of charged particles to the number of Dalitz pairs 
over all events of a given multiplicity class, and divide 
by the cocktail ratio to obtain the factor 1/e by which 
the raw data counts have to be multiplied. The discussion 
emphasises the importance to measure the reconstruction 
efficiency as a function of {Nch)- 



4.2.2 The Overlay-Monte-Carlo method 

A method of wider application is to reconstruct a MC- 
simulated sample of electron pairs from the sources un- 
der consideration, be it Dalitz decays with masses below 
200 MeV/c^ as discussed below, or the full hadronic cock- 
tail of electron pairs, always under condition of all accep- 
tance and analysis cuts and using standard analysis soft- 
ware. For realistic background conditions, the simulated 
pairs are overlaid one by one on data events of given {Nch) , 
and the desired correction factor is given by the average 
of the inverse probability of successful reconstruction, 




4.2.1 The Dalitz method 

Our sample of 'Dalitz pairs' with m <200 MeV/c^, but 
opening angles larger than 35 mrad and track pt above 
200 MeV/c is considered an excellent choice in place. It 
consists mainly of Tr'^-Dalitz pairs and a contribution of 
77-Dalitz pairs on the 10% level. The opening-angle cut 
reduces conversions to (5-10) % which can be corrected 
for. At present, we cannot exclude the (very interesting) 
possibility that a strong source of thermal electron pairs 
may compete on the level of a few percent j^. 

The required average inverse reconstruction efficiency 
is derived by dividing the measured number ratio of char- 
ged-particles to Dalitz pairs by the cocktail expectation,^^ 

With the event generator GENESIS 8¥ the number of 
electron pairs of the hadronic cocktail is calculated, but it 
is normalised to the number of neutral pions; using the 7r°- 
to-charged-particle ratio, we obtain the cocktail reference 
ratio in the required normalisation to charged particles, 

{nc^^fn^\M^ (4.2) 
{rich) yrich J {riT^o) 

For unit pair efficiency, all pairs that have both electrons 
in the acceptance of the spectrometer are reconstructed. 

By lower case n we denote numbers from physics simula- 
tion, by upper case A'^ measured counts. 



The pair efficiencies can be decomposed into the prod- 
ucts of the track efficiencies of the detectors for a mean- 
ingful check on the reliability of the simulation. Table 14.11 
compares full-track and pair efficiencies for a standard 
Dalitz sample with those obtained by piecewise multiplica- 
tion of individual detector efficiencies; all efficiencies were 
measured by the Overlay-Monte-Carlo technique. The close 

■^^ Since all components of the cocktail scale linearly with rich, 
ratios are constant. 





Reconstruction efficiency 
piecewise compound 


SiDCl 


94 % 


SiDC2 


88 % 


SiDC-tracks 


=^ 83 % 84% 


RICHl 


81 % 


RICH2 


86 % 


MWPC 


96 % 


full tracks 


=^ 55 % 55% 


pairs 


=^ 31 % 30% 



Table 4.1. Detector efficiencies separately and combined to 
full-track and pair efficiencies (left column) . Efficiencies of fully 
reconstructed tracks and pairs (right column). Overlay-Monte- 
Carlo technique, 1996 data sample. 
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Fig. 4.2. The multiplicity dependence of the efficiency to re- 
construct low-mass Dalitz pairs under acceptance and analysis 
cuts. The solid line is from Overlay-MC simulation. The data 
points show the decline in efficiency with increasing particle 
density assuming that the Dalitz yields are proportional to it. 
The integral efficiency was normalised to the MC value. Sta- 
tistical errors are comparable to or smaller than data circles, 
horizontal bars show binning of trigger distribution. 1996 data, 
Ref. EHI- 

agreement between the two indicates that the procedures 
of track reconstruction and pairing work properly. 

A further reduction in pair efhciency by about a fac- 
tor of two occurs in the course of the background rejection 
discussed in sect. 14.41 also these efficiency losses are un- 
derstood by the MC simulations. 

The Overlay-MC method applied to high-mass elec- 
tron pairs suffers from the low efficiency combined with 
acceptance losses which forbid to collect sufficient statis- 
tics samples. We have therefore based all corrections on 
products of track efficiencies constructed in such a way 
that on the one hand the correct dNch/drj dependence of 
the Dalitz pair efficiency as shown in Fig l4.2l is reproduced, 
and on the other, that the specific differences in pair ef- 
ficiency are taken care of in an approximate way. Those 
arise from the ^-dependent hit density and are treated 
separately. 

4.2.3 Nch dependent efficiency 

Reconstruction efficiencies for Dalitz pairs with mass be- 
low 200 MeV/c^ are displayed in Fig. 14.21 as a function of 
charged-particle rapidity density. The pairs are filtered for 
opening angles 0ee > 35 mrad to reduce conversions, and 
they obey the condition pf > 200 MeV/c on single-electron 
tracks. For the data points, the dependence on Nch is de- 
rived by postulating that sources scale with (Nch) while 
the absolute magnitude is fixed by normalising to the MC 
simulation as explained above. We observe a loss in MC 
pair efficiency by a factor of 2.2 when rapidity density is 
raised from 150 to 350. This degradation factor for the 
three '96-data analyses [551l56ll57| amounts to 2.30±0.13; 
the individual results lie remarkably close considering the 
large spread in absolute efficiencies. The causes of the ef- 
ficiency loss are twofold. For one, the recognition of hits 
in the SiDC's and of Cherenkov rings in the RICHes is 
deteriorated by an increasing number of background hits. 
In addition, with increasing track density also the chances 



increase that rejection of background tracks accidentally 
vetoes nearby signal tracks. 

The strong decrease of pair efficiency with {Nch) agrees 
very well with the Overlay-MC simulation also displayed 
in Fig. 14.21 In the '96 data analyses, this curve is used 
for efficiency correction on an event-by-event basis: signal 
pairs are stored with weights given by the inverse of the 
efficiency at given event multiplicity. The same procedure 
was applied in the '95 data analysis, but the (Nch) de- 
pendence is considerably weaker; the ratio quoted above 
is only 1.4 possibly because of a lower event back- 

ground compared to 1996, but also due to a different treat- 
ment of close hits in the SiDC's. 



4.2.4 (^-dependent efficiency 

With track density increasing towards small polar angles, 
the track reconstruction efficiency drops by more than a 
factor of three over the acceptance, as shown in the upper 
panel of Fig. 14.31 Since we will deal only with pair yields 
integrated over the acceptance, one is inclined to disregard 
this effect. 

A ^-dependent efficiency, however, combined with the 
limited 9 acceptance, affects pairs of different opening- 
angle characteristics in different ways; and we have seen 
already in Fig. 13. ll that 7r°-Dalitz pairs and pairs from uj 
decays have indeed very different 0cc behaviour. 

The acceptance condition affects wide open pairs con- 
siderably. For the given example of uj decays with typical 
0ee ~ 400 mrad, only pairs with both electrons at large 
6 are accepted and reconstructed with above-average ef- 
ficiency. The rise in pair efficiency for very large opening 
angles is seen in the lower panel of Fig. 14.31 The shallow 
minimum at Oee ~ 250 mrad is typical for the majority 
of the open pairs in our sample. We see also that the effi- 
ciency changes only by very little towards lower 0ce ■ The 
flat region even includes our sample of Dalitz pairs; owed 
to the standard cut 0ee > 35 mrad, their most-probable 
polar angle is as large as 190 mrad. So, compared to the 
'majority of open pairs' cited above, the reconstruction 
efficiency of the sample of Dalitz pairs is almost unde- 
graded. This is very satisfactory for the 'Dalitz method' 
described above. 

Yet, in order to protect characteristic features in spec- 
tra of pairs with opening angles comparable to the accep- 
tance of the spectrometer, the efficiency correction should 
involve the 9 distributions of the contributing tracks. We 
describe this method in the following. 

The 9 dependence in track efficiency was taken care 
of in an approximate way under the boundary condition 
to leave the measured efficiency for Dalitz pairs shown in 
Fig. 14.21 unchanged. The pair efficiency is factorised into 
track efficiencies which are modified according to 

£ (Nch, 9) = e.Zit (Nch) ■ f (9). (4.4) 

The function f{9) incorporates the measured 9 depen- 
dence of the track efficiency shown in Fig. 14. 31 upper part. 
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Fig. 4.3. Top: Track efficiency vs polar angle 9 drops as track 
density increases towards small 6. Bottom: Pair reconstruction 
efficiency vs pair opening angle 6>ee rises towards large ©ee 
where tracks of accepted pairs are forced to large 8. A shallow 
minimum occurs where tracks come closest to the minimum 9. 
Overlay-Monte-Carlo technique on '96 data. 

The procedure is to increment in an event of multiplic- 
ity Nch for each pair candidate a weight into the like-sign 
or unlike-sign spectrum which is 



' track 



1 



N any pair , j.j r. \ 
> track [^^ch,02) 



_ Dalitz 
-'pair 



{Nch) ■ f{ei) ■ f{02) 



(4.5) 



The function f{0) is suitably normalized to assure that 
the Dalitz sample is unchanged. 

The described event-by-event correction was not per- 
formed in the '95 data analyses while the '96-analysis of 
Ref. ISZ] was corrected for a truly two-dimensional effi- 
ciency. 



4.3 Pair acceptance 

The acceptance of the CERES spectrometer for electron 
pairs has been evaluated by MC simulations assuming 
uniform input distributions in mass m, pair transverse 
momentum pf^ , and rapidity y. How the CERES accep- 
tance depends on m and p'l'^ is shown in Fig. 14.41 for the 
standard rapidity range 2.1< r] <2.65; the normalisation 
is done cutting also the virtual photon distribution to 
this range in the input. We have applied here the stan- 
dard analysis cuts pt > 200 MeV/c and 0ce >35mrad. 
The results without cuts are shown for comparison. For 
masses below 100 MeV/c^, the acceptance without cuts 
approaches one, while the angular cut causes a steep re- 
duction. Higher up in mass, the single-electron pt cut re- 
duces the acceptance at low pair pj*^, acting like a mt cut. 
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Fig. 4.4. Acceptance of electron pairs under standard cut con- 
ditions (lower part) and without cuts (upper part) vs invariant 
mass, for five selected bins in transverse pair momentum pi" . 



while at pf > 400 MeV/c, the acceptance becomes more 
uniform and even, within a factor of two, independent of 
mass, due to some 'equalising' effect of the opening-angle 
cut. Still higher up, for m >500MeV/c^, the acceptance 
is essentially uniform both in mass and in p'^'^. 

An acceptance correction for pairs of 'new physics' 
sources is problematic in principle, since the decay dy- 
namics is not precisely known and may depend on addi- 
tional parameters, such as helicity angle. In comparison 
of theories to CERES data, proper accounts of the pair 
acceptance and other cuts are usually taken by the au- 
thors themselves who have the necessary insight into the 
dynamics of the modelled sources to do so. We therefore 
correct data only for reconstruction efficiency, not for pair 
acceptance. 



4.4 Optimising rejection of combinatorial background 

To foster confidence in the data analysis and to improve 
on data quality it was essential to understand the combi- 
natorial background in every detail. We have presented in 
the preceding section how the measured numbers of back- 
ground and signal pairs evolve with the rejection steps 
applied. Here, we use the tool of full MC simulation to 
arrive at a quantitative understanding of the background 
sources. It goes without saying that the elaborate sim- 
ulation of background rejection supplied a very effective 
handle to fine-tune rejection cuts. 

Figure E^l plots the number of electron tracks from var- 
ious physics sources which survive a given rejection step 
in the MC simulation. The decrease is either due to re- 
jection of single tracks that have unwanted properties like 
low Pt, bad matching or suspicious environment, or due to 
removal of recognised conversion or Dalitz pairs which are 
still intact and fulfil all conditions. 
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Fig. 4.5. MC simulation of number of electron tracks per Nch 
from various physics sources, along sequence of rejection steps 
and compared to measured background (full circles, 1996 data). 
Individual steps are 1: status after production (mild rejection, 
generous quality cuts), 2: tracking, 3: pt cut, 4: dE/dx cut, 5: 
fake rejection, 6: conversion rejection, 7: Dalitz rejection, 8: vr* 
rejection. 

The experimental points are numbers of background 
tracks from the post-production steps of the 1996 data 
analysis. They are derived from the measured total num- 
ber of background pairs accumulated in 42.2 10^ events, 
assuming Poissonian track statistics. The agreement be- 
tween absolute numbers from experiment and simulation 
is quite reassuring; also the size in relative suppression is 
in good agreement with the simulation. The residual back- 
ground remaining after all eight rejection steps is domi- 
nated by Dalitz pairs. 



5 Hadronic decay sources 

The 'conventional' sources contributing to the inclusive 
mass spectra of electron pairs in the mass range below 
1.5 GeV/c^ are free decays of light neutral mesons up to 
and including the (j). These contributions have been de- 
termined for p-Be and p-Au collisions with considerable 
precision in an experiment which combined the electron- 
pair spectrometer of CERES with the photon calorime- 
ter of TAPS • In addition to the inclusive electron- 
pair yield, the cross sections and pt distributions of the 
7r°,77, and u mesons have been measured via the electro- 
magnetic decay modes 77 and 7r°7, and an exclusive re- 
construction of the 7r° and ry-Dalitz decays 7r°, 77 ^ e+e^7. 
These data, supplemented by older measurements from 



NA27 IHSj, allowed us to simulate the electron-pair mass 
spectra originating from decays of neutral mesons. The re- 
sult was subject to all experimental cuts and folded with 
mass resolution. Within error limits it was concluded that 
the inclusive mass spectra of electron-pairs measured in 
450 GeV/c p-Be and p-Au colhsions are consistently de- 
scribed by the expected e~^e~ pair 'cocktail' of neutral 
meson decays IT6| . 

The hadronic cocktail serves as our reference standard 
of hadronic sources to expose effects which are specific to 
nucleus-nucleus collisions, i.e. spectral shapes and yields 
which cannot be described by a mere superposition of 
nucleon-nucleon interactions. Supported by evidence from 
light collision systems 86 , particle ratios were originally 
assumed constant from p-p to heavy systems [ST], and 
particle yields were scaled with charge multiplicity. This 
cocktail was instrumental in gauging the low-mass excess 
in 200 GeV/c S-Au collisions QJj and in 158 GeV/n Pb- 
Au collisions from the first (1995) data set 

The hadronic cocktail is calculated with the Monte 
Carlo event generator GENESIS [HJ described in consid- 
erable detail elsewhere [T^. It has been improved since 
then but the changes with respect to previous CERES 
analyses of Refs. |48ll49ll5()ll51j are subtle. These publica- 
tions are superseded by the present paper. 

We sketch here only the main points. Mesons are pro- 
duced with cross sections scaled up from p-p or p-A, tak- 
ing into account the modified pt and rapidity distribu- 
tions. Open charm contributions have been neglected on 
the basis of reliable estimates jS^- Mesons are then al- 
lowed to decay with known branching ratios. Decay kine- 
matics are simulated for the Dalitz decays PS e^e~^ 
of the pseudo-scalar mesons 7r°,77 and rj' , for the direct 
decays of the light vector mesons, V e+e~, and for 
the a; Dalitz decay lo — > e^e~-K° . Electron momenta are 
Lorentz-transformed into the laboratory system and con- 
voluted with the experimental resolution profile. The sim- 
ulated events are subject to the same filters concerning 
acceptance, pt and opening angle. 

Production cross sections of the neutral mesons, their 
rapidity and transverse-momentum distributions are es- 
sential ingredients for simulating the decay contributions 
to the dilepton spectrum. When data on hadron produc- 
tion from Pb-beam experiments at the SPS (NA44, NA49, 
NA50, WA98) became available in time for the 1996 data 
analysis, we used this input from Pb-beam data whenever 
possible [Knil51ll56| . Information not directly available was 
derived from the statistical model which describes ratios of 
integrated hadron yields at chemical freeze-out very well 
with only two fit parameters, the temperature and the 
baryon chemical potential; the particular values used are 
T= 170 McV and = 68 MeV 89 . The coUision system 
is modelled in a state of collective transverse expansion 
which is based on the observation [501 that inverse-slope 
parameters T of transverse momentum spectra, except for 
pions, systematically increase with mass '91'. 

Data on pt distributions of pions have been measured 
by CERES NA49 NA44 ^5], and WA98 

They are exceptional in that they can not be described by 



26 



G. Agakichiev et al.: e''"e -pair production in Pb-Au collisions at 158 GeV/nucleon 



decay 


crtot/o-„o^tot 


branching- 
ratio 


To [MeV] 


Form factor 


tt" e^e~7 


1. 


(1.198 ±0.032) X 10"^ 


100 / 230 


pole approximation 
&=5.5 GeV"^ 


■q — -> e'^e^'j 


0.085 


(6.0 ±0.8) X 10"'' 


238 


pole approximation 
6=(1.9± 0.4) GeV'^ 


p e^e~ 


0.094 


(4.67 ±0.09) X 10~" 


263 




Lu 7r°e+e~ 


0.069 


(7.14 ± 0.13) X 10"'= 
(5.9 ± 1.9) X 10"^ 


265 


Breit-Wigner 
mp=0.65, rp=0.05 GeV 


(f) — > e+e~ 


0.018 


(2.98 ±0.04) X 10"" 


292 






0.0078 


^ 5.6 X 10"* 


285 


Breit-Wigner 
mp=0.76, rp= 0.01 GeV 



Table 5.1. Summary of the parameters used in the generator to simulate the contributions from hadron decays. 



a single exponential. The tt" spectrum is generated with 
two slopes: at mt < 200 MeV where it is dominated by sec- 
ondaries, with To— 100 MeV, and for the higher part with 
an inverse slope of To— 230 MeV. The inclusive mt dis- 
tribution for neutral pions measured by WA98 is extrap- 
olated to small mt using the charged pion distributions 
from NA44. At small mt the spectra are dominated by sec- 
ondary decays from heavier mesons. The decay r] —f 3n° 
is added separately by hand to the n° distribution. 

The rapidity distribution |97| of negatively charged 
hadrons, described as a Gaussian centred at ymax= 2.9 
with ay— 1.5, has been adopted for all mesons. While the 
widths of the hadron rapidity distributions decrease with 
particle mass in proton induced collisions, this is not ob- 
served in lead-induced collisions [55].^^ 

The parameters from Pb-beam data and the statistical 
model are given in Table IsTI Meson production ratios are 
implemented relative to the number of 7r°'s and include 
feeding from heavier resonances. Compared to the earlier 
reference to proton induced coUisions |l(i| . the production 
ratios of heavier hadrons are enhanced as, e.g. the ti/tt^ 



j99 and the 20/ (tt^ 



ratio 100,101 . The meson 



yields are normalised to the charged-particle density by 
the ratio 

{N^o)/{N,h) = 0.44. (5.1) 

Brackets denote averaging over the CERES acceptance. 

The decay branching ratios given in Table IFTI are from 
Ref. (102) . To simulate the Dalitz decays, the KroU-Wada 
expression jl03| is multipHed by the electro-magnetic tran- 
sition form factors fitted to the LEPTON-G data [TUl] . 
The pole approximation F{M'^) = (1 — bM'^)~^ is used 
for the determination of the form factors of the 7r° and 
the r] Dalitz decays |l()4j . For the uj and the ry', the form 
factors are determined by fitting a Breit-Wigner function 



|F(M2)|2 



(M2 - m2)2 + mir^ 



(5.2) 



to describe the resonant behaviour according to the vec- 
tor dominance model. The parameters used are listed in 



Particle ratios taken at mid-rapidity are therefore the same 
as those from partially or fully integrated yields. 



Table l^m The direct decays of the vector-meson were gen- 
erated following Gunaris and Sakurai |105j . 

The 2-body decay of the p meson has been re-evalua- 
ted HHI for the new GENESIS code HHZl- We defer details 
of the revised formula as it is implemented in the 2003 ver- 
sion of the New GENESIS to the appendix. The resulting 
mass distribution, due to a Boltzmann-type phase space 
factor e~*^/^ and a momentum dependent phase space, 
both omitted in the previous code, receives a shoulder on 
the low-mass side and a steeper falloff to higher masses. 

All decays where assumed isotropic in the rest frame of 
the decaying meson except for the Dalitz decays to e+e~7 
which follow a 1 ± cos^(6') distribution, where 9 is mea- 
sured with respect to the virtual photon direction. 

Figure IHTI depicts the cocktail based on the Pb-beam 
data with all corrections. It is the mass spectrum of pairs 
in the CERES acceptance 2.1 < 77 < 2.65 with standard 
cuts Oee > 35 mrad and pt > 200 MeV/c . The generator 
output has been folded with the mass resolution function 
of the 1996 data set, which includes both the momen- 
tum resolution as the dominant source of smearing as in 
eqn. (3.9), as well as the resolution in pair opening angle; 
the latter is approximately 



(72, 



ere) 



(5.3) 



where ag ~ 0.6 mrad and ~ 3.0 mrad are the angu- 
lar track resolutions taking information from all detectors 
together. The effects of bremsstrahlung emission by elec- 
trons traversing the detector material have been included 
in the simulations, but due to the low material budget in 
the acceptance of X/Xo « 1 %, they are hardly notice- 
able. The overall mass resolution is about 6% in the p/w 
region and 7% in the region of the (see sect. 3.3). 

The cocktail for Pb-Au coUisions is compared in Fig. 15. II 
to the cocktail of hadronic sources up-scaled from p-p in- 
teractions; the latter is also folded with the CERES reso- 
lution obtained for the '96 data analysis. 

The systematic errors of the normalised decay cocktail, 
relevant for the numerical comparison to the experimental 
pair data in sect. 7, are discussed separately for the mass 
regions below and above 200MeV/c^. 

The low-mass region is dominated by the 7r°-Dalitz 
decay, with 7]-Dalitz contributing about 15%. Errors arise 



G. Agakichiev et al.: e^e -pair production in Pb-Au collisions at 158 GeV/nucleon 27 



p-p particle ratios 
Pb-Pb particle ratios 




of 450-750 MeV/c^, and 6% for m >750MeV/c2. Taking 
these supposedly independent sources of uncertainty to- 
gether, an overall systematic error of 30% is estimated as 
the weighted average for the high-mass region. 

The integral number of electron pairs per charged par- 
ticle in the CERES acceptance expected from hadronic 
sources under standard cut conditions > 200 MeV/c 
and 0ee > 35 mrad is 



(dNee/dm) 



Genesis 



(9.27 ±0.74) 10-^ 
(2.27 ±0.7) 10-^ 



m < 
m > 



200 MeV/c2 
200 MeV/c^, 



(5.4) 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

m,, (GeV/c') 



Fig. 5.1. New GENESIS hadronic cocktail based on Pb-beam 
data and the thermal model compared to GENESIS based on 
p-p collisions. Shown are mass spectra from simulated neutral 
mesons decays with acceptance, pt and opening-angle cuts and 
after folding with 1996 mass resolution. 



from the relative production cross section of the it° and 
from the parametrisations of the input rapidity and trans- 
verse momentum distributions. Normalising the decay cock- 
tail to the number of charged particles is a very pow- 
erful means to keep the error in the ratio (71-^0 /rich) of 
eqn. (4.26) small; it is estimated to be at most 5%. The 
rapidity distribution is uncritical since it is taken directly 
from pion data; the errors should be less than 3%. The 
transverse momentum distribution, due to the single-electron 
Pt cut, is a little more critical. However, quite different as- 
sumptions on the shape of the pt spectrum above the cut, 
e.g., using from CERES rather than tt^ from NA49, 
lead to differences in yield of only very few percent. This 
error is therefore also estimated to be below 5%. Assum- 
ing that the systematic errors are uncorrelated, a total 
systematic error of 8% is obtained for the low-mass re- 
gion. 

In the high-mass region, the errors are dominated by 
those of the relative production cross sections of the higher- 
mass mesons and of the detailed properties of the electro- 
magnetic decays. Since most of the particle yields rele- 
vant for the cocktail have not directly been measured, or 
like the 0, suffer from experimental controversy 100,101 , 
the statistical model predictions have been used instead, 
and it is therefore their uncertainties which enter. Judg- 
ing the average fit quality to measured particle ratios 89 , 
we estimate these errors to be 20%. The uncertainties in 
the branching ratios and in the transition form factors 
have been discussed in detail in Ref. ^Hl- They contribute 
«15% for m <450MeV/c^, «30% in the mass range 



where we have quoted the systematic errors of 8% and 
30% for the two mass regions, respectively. 



6 Results 

6.1 Samples of reconstructed pairs 

For an overview, the pair samples reconstructed from the 
two data sets of 158 GeV/n Pb-Au collisions taken by the 
CERES Collaboration in 1995 and 1996 are listed in Ta- 
ble 16.11 All analyses apply standard cuts on track pt and 
and pair opening angle 0ee- 

The table hsts the pair signal (S) obtained by subtract- 
ing the combinatorial background (B) from the measured 
numbers of unlike-sign pairs, for the mass ranges below 
(Dalitz pairs) and above 200 MeV/c^ (open pairs). In this 
Table, pair yields are not corrected for reconstruction ef- 
ficiency (other than in the following figures). In analysis 
no. 5, event mixing was employed to obtain combinatorial 
background of unlike-sign pairs. Quoted errors are abso- 
lute statistical errors in the respective sample numbers 

amounting to as = \J + B. The pair reconstruction 

efficiencies are given in the last column. 

We note a considerable spread among the results of 
different analysis efforts. This originates from the different 
values of the rejection cuts used along the various analysis 
chains. However, normalised pair yields, i.e. numbers of 
pairs per event, per charged particle, and corrected for 
pair efficiency, as will be shown below, are very stable: 
the relative spread in the number of Dalitz pairs turns 
out to be less than 15%. 



6.2 Inclusive mass spectra 

Figure 16.11 shows the mass spectrum of Ref. |49II54| from 
1995 together with that of Ref. IKUKTll^ from the 1996 
data set. The trigger ccntralities correspond to the most 
central 33% and 26% of the inelastic cross section for 
1995 and 1996, respectively. The signal is obtained by sub- 
tracting the smoothed like-sign pair background from the 
spectrum of unlike-sign pairs. The differential pair yield 
(dNee/dmee) psr event is normalised to the mean number 
{Nch) of charged particles in the acceptance. The brack- 
ets indicate particle yields per event measured within the 
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Work 


Events 


Dalitz 
Pairs 


S/B 


Open 
Pairs 


S/B 


Open/Dalitz 
Pair Efficiency 


Year of run 


{dN,h/d?7) 












1. 1995a 
Refs. 49,54 


8.55-10" 
220 


1038±55 


1.05 


648±105 


1/8 


10.9/11.0% 


2. 1995b 
Ref. 55 


8.55-10" 
220 


1044±53 


1.0 


468±104 


1/11 


8.0/10.9% 


3. 1996a 
Refs. 50,51,56 


42.0-10" 
250 


5631±129 


1.03 


2018±237 


1/13.4 


8.8/11.9% 


4. 1996b 
Ref. 55 


42.0-10" 
250 


5013±128 


0.88 


1722±236 


1/15.7 


9.5/10.8% 


5. 1996c 
Ref. 57 


41-10" 
250 


3537±103 


1.04 


1305±194 


1/13.9 


4.3/ 5.6% 



Table 6.1. Analysis results of the 158 GeV/n Pb+Au data sets of 1995 and 1996 with total number of events analysed and 
mean charged density of trigger. Listed are the signals S ± as of reconstructed 'Dalitz' and 'Open' pairs of mass below and 
above 200 MeV/c^, respectively, after background subtraction. Signal-to-background ratios S/B are quoted. Pair reconstruction 
efficiencies (l/e)^^ from open-pair and Dalitz-pair analyses are given in the last column. Yields are not efficiency corrected. 
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Fig. 6.1. Inclusive e^e~ mass spectra from '95 and '96 data 
analyses. The pair yield is efficiency corrected and normalised 
to the number of charged particles in the acceptance. Standard 
analysis cuts on track pt and pair opening angle have been 
applied. Data correspond to the top 33% and 26% of ainei 
for the '95 and '96 data sets, respectively. Vertical bars give 
statistically independent bin-to-bin errors. 



CERES acceptance. Reconstruction efficiency has been 
corrected event- by- event by weighting pairs with the value 
of the inverse pair efficiency at the particular centrahty. 
Shown are the bin-to-bin statistical errors of eqn. (3.18). 
To include the normahzation errors discussed in sect. 13.71 
as it may be relevant when individual data points are com- 
pared to other data, as in this figure, or to theoretical- 



model predictions, the errors should be multiplied by 1.4. 
Note, however, that these larger errors are no longer sta- 
tistically independent. 

It is apparent that the '95 data points lie systemati- 
cally higher than the '96 data points. Apart from the small 
difference in mean trigger centrality there is no relevant 
change in setup. The '95 analysis |54II48II49| followed the 
strategy to optimise the quantity Seff ~ / B. But while 
the signal efficiency e was determined by Monte-Carlo, the 
background B was taken as the measured like-sign sample; 
this strategy was discarded when the full MC simulation 
became available for the '96 analysis. Although we could 
find absolutely no indication that statistical ffuctuations 
in B might have steered the analysis towards a 'better' 
final sample, such possibility cannot be strictly excluded. 

The re-analysis of the '95 data ^55^ resulted in a mass 
spectrum closer in absolute yield and shape to the '96 
mass spectrum shown in Fig. 16.11 In this analysis, any 
involuntary bias was avoided by sampling the distribu- 
tions in a random automatic variation of all cut settings 
simultaneously, and then choosing the centres of gravity. 
However, seeing no direct evidence for a biased tuning of 
cuts, we have chosen to keep the original '95 analysis to 
expose our actual systematic uncertainties. We return to 
this issue below. 



6.3 Centrality dependence 

The efficiency-corrected yield of Dalitz pairs is a solid ref- 
erence for linear N^h dependence. We have a first look 
at the centrality dependence of the open pair yield by 
comparing in Fig. 16.21 the centres of gravity of three N^h 
distributions: that of triggered events without any further 
condition, and two others taken from events which contain 
a Dalitz or an open pair candidate, respectively. Because 
of the large open-pair background, the signal is obtained 
by subtracting the Nch distribution of like-sign-pair events 
from that of unlike-sign-pair events. The A'c/i-dependent 
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Fig. 6.2. The dNch/drj distribution of the trigger (full circles), 
and of events containing a Dalitz (open circles) or an open pair 
candidate (filled squares). Centres of gravity are marked by the 
respective symbols in a box at the bottom. 



pair detection efficiency has been corrected for using the 
curve of Fig. 14.21 

We see from Fig. 16.21 that the centres of gravity of 
Dalitz and open-pair samples are shifted progressively up- 
ward; the mean multiplicities are 285 and 310, respec- 
tively, compared to {Nch)=250 of the trigger distribu- 
tion. The averages of the trigger distribution calculated 
for events depending linearly or quadratically on Nch are 
287 and 315, respectively. We conclude that the increase 
of open-pair production is better described by a quadratic 
dependence on Nchi than by a linear dependence. We will 
have a closer look at the differential centrality dependence 
in sect.Elbut note that the current finding bears no refer- 
ence to the hadronic cocktail. 



6.4 Invariant transverse-momentum spectra 

Considerable physics potential resides in the spectra of 
invariant transverse pair momentum, pf'^, as displayed in 
Fig. 16.31 for the two data samples and three mass bins. 
The spectra are normalised and corrected for reconstruc- 
tion efficiency. In the Dalitz region, at low the spec- 
trum is void due to the analysis cut; both samples agree 
very well. Data points scatter considerably for the two 
other mass bins. The full physics relevance of these spec- 
tra emerges only in comparison to the expectations for 
hadronic sources. 



6.5 Stability of results 

The optimisation of quality and rejection cuts was done 
with no feedback from the signal itself not to exploit sta- 
tistical fluctuations likely in a sample that small. Still, 
at a signal-to-background level of order 1/10 one is con- 
cerned how reliable and stable a signal can be that results 
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Fig. 6.3. Invariant pair transverse momentum spectra for 
three mass ranges and both data sets ('95 squares, '96 circles). 
Statistical errors only. Trigger centrality 33% and 26% for the 
'95 and '96 data sets, respectively. 



from subtracting two almost equal large numbers. Even 
by careful inspection of raw mass spectra by shape, or by 
magnitude, it would be quite hard to distinguish above 
200 MeV/c^ the signal from background. 

Because of cuts and kinematics, it so happens that the 
combinatorial background spectrum peaks around m « 
2p('"*/c= 400 MeV/c^ and is surprisingly similar in shape 
to the low-mass enhancement. 

This raises the question, whether combinatorial back- 
ground has been insufficiently subtracted. Aside from care- 
ful studies of analysis cuts, rejection cuts, the determina- 
tion of the reconstruction efficiency as a function of cen- 
trality, we like to address this question in a quantitative 
way. A most likely cause leading to wrong subtraction of 
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Fig. 6.4. Results of all CERES electron-pair analyses for 
158 GeV/n Pb-fAu. Plotted are normalised pair yields inte- 
grated from 200 MeV/c^ until 1.1 GeY/c^ vs background-to- 
signal ratio. Errors are statistical. Slashed entries refer to re- 
laxed rejection or quality cuts in otherwise identical analyses. 
Enclosed numbers refer to Table l(T!n 
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background is an undetected asymmetry in reconstruc- 
tion efficiencies, i.e. for like-sign as compared to unlike- 
sign pairs. It is shown in Appendix C that a 5% asym- 
metry is required to fake an apparent enhancement factor 
of three over the hadronic sources, assuming the S/B= 
1/13 situation of the '96 analysis. However, this asymme- 
try was measured to be less than 1%, with confidence limit 
of 90% as reported in sect. 13.61 

If the low-mass enhancement is faked by leakage of 
some amount of combinatorial background into the spec- 
trum of signal pairs, the enhancement should increase with 
the amount added. Figure 16.41 summarising all CERES 
pair analyses for 158 GeV/n Pb-Au, does not show that. 
Rather, we see a pair signal which is, in view of the er- 
rors, surprisingly stable despite very large variations in 
background level. 



of 18%'^'' is larger than the rms deviation of Dalitz pair 
yields from all our analyses which amounts to 12%. 
The normalised yields of signal or open pairs 

.1.1 

AA('Open') = / dm (dNee(m)/dm)/(Neh) (6.2) 

JO.2 

plotted in Fig. 16. 41 displav a relative spread of 24%, twice 
that of the Dalitz sample and now well above the sys- 
tematic Dalitz-sample error. We presume that the larger 
systematic uncertainties in open pair yields are caused 
by the large combinatorial background. We take the rel- 
ative sample error of 24% in number of normalised open 
pairs as a reliable measure of the magnitude of systematic 
uncertainties, arguing that they were obtained in rather 
independent efforts and by using diverse strategies. 



6.6 Systematic errors 

Estimates of systematic errors in the mass-integrated, nor- 
malised yield of 'Dalitz' pairs^^, 

AA('Dalitz') = / dm (dNee(m)/dm>/(Nch), (6.1) 
Jo 

is given in Table The definition of the 'Dalitz' sam- 
ple includes the opening angle cut of 35 mrad and the pt 
cut at 200 MeV/c. The largest contributions of about 10% 
each arise from detector efficiencies and the triggered cen- 
trality (Nch), while the uncertainties due to matching and 
rejection cuts account for about 7%. The resulting total 



Note that the use of calligraphic Af for normalised pair 
yields is restricted to this section. 



7 Final Results and Comparison to Hadronic 
Sources 

7.1 Inclusive mass spectra 

The invariant mass spectrum of e^e~ pairs produced in 
158 GeV/n Pb-Au collisions as obtained from a weighted 
average of the individual data sets, henceforth called 'uni- 
fied mass spectrum', is shown in Fig. 17.11 The squares 
of the inverse relative statistical errors have been used 
as weights. Plotted is the differential yield per event cor- 
rected for pair efficiency and normalised to the number of 
charged particles in the acceptance, (Nch) ~ 133, which 
corresponds to the rapidity density {dNch/drf) « 245 and 
the most central 28% of the reaction cross section. 

If errors from both tracks would be correlated, they would 
add to a total of 24% instead. 



Source of systematic error 


Comment 


Relat. error 
on Noc/Nch 


I. 

MC efficiency correction 
SiDC's and Pad Chamber 
RICH'S 


3.0 % per hit on track 
5.0 % per ring on track 


± 7.3 % 
± 10.0 % 


Sum detector efficiencies 




± 12.4 % 


II. 

Efficiency of analysis cuts 
11 matching cuts 
2 quality cuts 
7 rejection cuts 


1 % per cut on track 

2 % per cut on track 
1 % per cut on track 


± 4.7 % 
± 4.0 % 
± 3.7 % 


Sum analysis cuts 




± 7.2 % 


III. 

Triggered mean Nch 
No of events analysed 


comparison to other 
SPS experiments 


± 10.0 % 
± 3.0 % 


Total 




± 17.7 % 



Table 6.2. Estimates of systematic errors in normalised yields of Dalitz pairs. The uncertainties in efliciency for the two tracks 
of the pair are assumed to be ancorrelated and are added in quadrature like the contributions from different cuts or detectors. 
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Fig. 7.1. Unified mass spectrum corresponding to average 
trigger centrality of 28% cr/cTinei- Yields are corrected for re- 
construction efficiency and normalised to (Nch)- Statistical bin- 
to-bin errors are shown by vertical bars, systematic errors by 
brackets. Horizontal bars indicate size of mass bins. Also shown 
is cocktail of hadron decays and its sum (solid line) . 



The mass spectrum is compared to the expectation 
of hadron decays modelled with the updated GENESIS 
code. The resonance structure of the light vector mesons 
is hardly visible in the measured spectrum, and within 
statistical errors, the data points are compatible with a 
smooth curve. Note that the resonances were clearly visi- 
ble in the p-Be and p- Au invariant mass spectra(Fig. 
despite a comparatively poor resolution of 9% at p/cu. 
Note also that the p/uj resonance seems to become vis- 
ible above the continuum in the mass spectrum selected 
for Pf"" > 500 MeV/c^ of Fig. |721 to be discussed below. 

In the (p region, the expected mass resolution of 7% be- 
ing potentially sufficient, statistics hampers further con- 
clusions. We nevertheless quote the ratio of the observed 
yield (in the mass region 0.95< m <l.lGeV/c^) rela- 
tive to the cocktail which is dominated by the statistical- 
model result for the (j) in this region. We obtain a value 
of 0.68±0.28. With one standard deviation below unity, 
the statistical accuracy is insufficient to settle the pending 
controversy on the (j) |1(JUII101| . 



7.2 Enhancement factors 

All CERES results for Pb-Au collisions show pair yields 
in the 7r°-Dalitz region which are in good agreement with 
predictions from known hadron decays. For masses above 
200 MeV/c^, however, the data overshoot the expectation 
from hadron decays significantly. The largest enhancement 
over the hadronic cocktail is in the mass range between 400 
and 600 MeV/c^ where it reaches a magnitude of six. 



Integrating the normalised yields up to 200 MeV/c^, 
for the fraction dominated by 7r°-Dalitz pairs (A), and 
from 200MeV/c^ upward for open pairs (B), we obtain 



(8.52 ± 0.20 [stat.] ± 1.54 [syst.]) x 10"^ (A) 
(5.25 ±0.43 [stat.] ±1.26 [syst.]) x 10"^ (B) ' 



(7.1) 



We have quoted the statistical and the systematic errors 
of 18% and 24%, respectively. 

The ratios of the measured data to the integrated yields 
of the decay cocktail (given by eqn. (5.33)) for the two 
mass regions are 



{N,,)/{N,h) 



decays 

0.92 ± 0.02 [stat] ± 0.17 [syst] ± 0.07 [decays] {A) 
2.31 ± 0.19 [stat] ± 0.55 [syst] ± 0.69 [decays] \b) 



where our estimate of the systematic error in the decay 
cocktail is given separately. For the low-mass continuum 
region 200 < m < 600 MeV/c^, the enhancement is even 
larger, 

J^ = 2.73±0.25 [stat] ±0.65 [syst] ±0.82 [decays]. (7.3) 



7.3 Centrality dependence 

The enhancement factors J-' are plotted vs charged-particle 
density for three mass regions in Fig. 17.21 The enhance- 
ment for the mass region 200< m < 600 MeV/c^ reaches 
about 4 at the most central collisions, and it is seen to rise 
about linearly with charged-particle density; a straight- 
line fit to the five data points gives a slope value deviating 
from zero by 4 standard deviations. This establishes that 
the integrated pair yield itself has a stronger-than-linear 
dependence on charged particle density. (By construction, 
the cocktail yield per Nch does not depend on Nch)- This 
enhancement extends to the resonance region, although 
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Fig. 7.2. Enhancement factor vs dNch/dy for three mass 
regions, unified data. The horizontal line indicates 1 of 
hadronic sources. The centre plot exhibits the largest effect, 
and the data points are consistent with a straight line that 
passes at dNch/drj = through '1' of the ordinate. 
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Fig. 7.3. Unified mass spectra for selected contiguous ranges 
in Nch, {Nch)= 150, 210, 275, and 370 (from upper left to lower 
right). Unified data, standard analysis cuts. 



with reduced significance and smaller values of the en- 
hancement factor. Unfortunately, statistical uncertainties, 
do not allow to confirm or refute interesting details in the 
growth of the normalized yield with (Nch), like a threshold 
effect or a saturation behaviour. 

To see how the enhancement-typical spectral shape 
evolves with increasing centrality, unified mass spectra 
for different ranges in Nch are shown in Fig. 17.31 The se- 
ries of spectra clearly demonstrates that the excess over 
the cocktail for all four spectra occurs between the two- 
pion threshold and the p/w-resonance position, increasing 
with centrality. The spectra, being statistically indepen- 
dent samples, corroborate that the largest enhancement 
is around 500 MeV/c^, well below the p/ui position. With 
respect to the central issue of in-medium modifications, 
one should be careful not to establish a direct link to the 
enhancement factor as we defined it. The reason is that 
pion annihilation per se, i.e. with a vacuum p, proceeding 
in the hot fireball is already a large, if not the dominant, 
contribution to the enhancement factor. 



7.4 Pair transverse momentum 

The Lorentz-invariant transverse momentum distributions 
of produced electron pairs, observed first in the 1995 data 
analysis are shown in Fig. l7.4l for the combined 95/96 




Fig. 7.4. Unified invariant pair-pt spectra me^O.2 (left), 0.2< 
me(< 0.6 (centre), and ruee > 600MeV/c (right). 28% (j/ainei- 



data for the 7r°-Dalitz region and for open pairs. At low 
pair pt, denoted by pf'^, the yield is suppressed due to the 
single-track pt cut at 200MeV/c. As to be expected, there 
is good agreement between data and hadronic cocktail in 
the Dalitz region. However, the information from the pair 
Pt distribution in the open-pair mass range is striking: 
the enhancement grows towards small pair pf'^ despite the 
Pt >200 MeV/c condition on single electron tracks. The 
surplus originates from decays of virtual photons within 
the fireball that are created favourably at rest. The en- 
hancement above 500 MeV/c pair transverse momentum 
is considerably reduced. Whether this observation trivially 
reflects the annihilation kinematics in a thermal pion gas, 
or contains information about modified hadron properties, 
will be discussed in sect. 8. 

Conversely, the size of the average transverse pair mo- 
mentum has a remarkable influence on the shape of the 
inclusive mass spectra. We observe in Fig. 17.51 that the 
measured yield dramatically overshoots the hadron decay 
contributions when low pair momenta are selected. The 
enhancement in the mass range 500-700 MeV/c^ reaches 
locally an order of magnitude. In contrast, for larger trans- 
verse pair momenta, the measured yield comes pretty close 
in shape and magnitude to the cocktail expectations; the 
resonance region becomes visible as it gains in yield, even 
on absolute scale, while the previously amplified contin- 
uum at 400< TO <600 MeV/c^ is being deflated. 



8 Physics discussion 

The two CERES runs of 1995 and 1996, combined to 
the largest statistics sample of 158 GeV/n Pb-Au colli- 
sions taken so far, corroborate the enhancement in low- 
mass electron-pair production over that from hadron de- 
cays originally observed in 200GeV/n S-Au collisions by 
CERES 1:2: and in S-W collisions by Helios-3 U. Two 
additional CERES Pb-Au runs have been performed since 
then with the new TPC and improved mass resolution, 
one at the reduced energy of 40 GeV/n [Knill()8lll()9| the 
results of which are shortly addressed below, the other at 
158 GeV/n [ITn) with quantitatively consistent results. 



G. Agakichiev et al.: e^e -pair production in Pb-Au collisions at 158 GeV/nucleon 33 




Fig. 7.5. Unified mass spectra for two selected ranges of transverse pair momentum, pf*^ < 500 MeV/c (left), and 
pT > 500 MeV/c (right). Centrality 28% a/ainei- 



In the region of the 7r°-Dalitz decay, m < 200 MeV/c^, 
the yield of the unified sample agrees with the cocktail 
within the estimated systematic uncertainties. In the mass 
range upward of 200 MeV/c^, the normalised pair yield 
of the unified data analysis significantly exceeds the yield 
from hadron decays per charged particle, and the enhance- 
ment factor is given together with our estimates of statis- 
tical and systematic errors in sect. 17.21 

The systematic errors of the decay cocktail were dis- 
cussed in sect. 4. To address the reliability of the generator 
in physics terms, the situation has remarkably improved 
since the practice of scaling yields up from p-p to Pb-Au 
collisions was abandoned in favour of using particle ratios 
from statistical model systematics based on Pb-beam data 
itself. Still, we shortly recall here some particle ratios that 
have become subject to speculations of being grossly en- 
hanced (i.e. by factors more than two) in nucleus-nucleus 
compared to p-p reactions, even beyond the statistical- 
model systematics. 

First, the rj/ir" ratio received interest from the fact 
that the ry-Dalitz decay is the most important single com- 
ponent of the hadronic cocktail being suspect of enhanced 
production. It was argued however, that anrj/n" ra- 

tio sufficiently large to explain the electron pair enhance- 
ment had not been observed in the photon yield measured 
by WA80 [112 . 

In the region below the p/to, the gap between cock- 
tail and data might be filled by raising the cj-Dalitz de- 
cay contribution as has been proposed by V. Koch some 
time ago We hold against that there is no indica- 

tion otherwise for a strongly enhanced to production. In 
addition, a boost in w-Dalitz production is limited to an 
upper margin set by the data of Fig. IT.ll to the direct de- 
cay Lu e+e^. Eventually, there is absolutely no reason 
to believe that open-charm production in nucleus-nucleus 



collisions should be as enormously enhanced as to explain 
the low-mass enhancement jSHj- 

We return to state that the enhanced production of 
low-mass electron pairs cannot be attributed to decays of 
produced hadrons. The excess has to originate from pro- 
cesses which are active during the lifetime of the fireball, 
i.e. between the onset of hadronisation and kinetic freeze- 
out, if of hadronic origin, and before hadronisation, if cre- 
ated during the plasma phase. Present theoretical studies 
allocate only a very minor fraction to the partonic part 
which reflects the supposedly small 4-volumes of decon- 
fined matter at SPS energies. 

After having tried to give an overview of theoretical 
models of dilepton production in the introduction, this dis- 
cussion will be guided by a few representative theoretical 
models; the spectral function approach 6 , the dropping- 
mass scenario T , and pion annihilation with a vacuum p. 
In all calculations to be shown, the same fireball model 
has been used to describe the space-time evolution |114) . 

Which are the experimental signatures of pion annihi- 
lation, the process most widely ascribed to take place in 
the dense hadronic fireball ? As for any binary process, the 
annihilation rate is expected to scale with the squared den- 
sity of the particles annihilating in the fireball (of course, 
the argument applies as well to qq annihilation). Unfor- 
tunately, the centrality dependence of the dilepton yield 
is a topic barely addressed by full transport calculations, 
except jnSI- 

We take the observed stronger-than-linear scaling of 
the pair yield (in the mass region of the strongest en- 
hancement) with Nch as strong evidence in support of the 
two-body annihilation reaction. This deserves some words 
of justification. An increase in Nch signals a larger pion 
density only to the extent that it is not compensated by 
an associated increase in volume, maybe even in lifetime. 
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Fig. 8.1. The pt" spectra of Fig. 17.41 compared to (i) free hadron decays without p (thin solid line), (ii) model calculations 
with a vacuum p spectral function (thick dashed line), (iii) with dropping in-medium p mass (thick dashed-dotted line), (iv) 
with a medium- modified p spectral function (thick solid line). For a cocktail including the p decay see Fig. 17.41 



of the fireball. Therefore, pion annihilation does not neces- 
sarily go along with quadratic scaling in Nch- Conversely, 
however, an observed stronger-than-linear scaling of the 
pair yield with Nch is sufficient reason to infer a binary 
reaction at work, i.e. strongly suggesting pion annihila- 
tion in hadronic matter, or qq annihilation in the plasma 
phase. 

As other annihilation processes in a thermal medium, 
TTTT annihilation takes place favourably at small relative 
momentum of the constituents. This behaviour was in- 
deed encountered already in the pair transverse momen- 
tum spectra of Fig l7.4l From Fig l8.1l shown here, we learn 
in addition that pion annihilation with a vacuum p is hard 
to distinguish, by yield and shape of its p^"^ spectrum, from 
the medium-modified spectral function and the dropping- 
mass approaches. 

The invariant mass spectra tell more about the physics 
processes involved. In Fig. 18.21 the measured invariant 
mass spectrum is compared to the model calculations. It is 
evident at first sight that pion annihilation with a vacuum 
p (thick dashed line) does not describe the shape of the 
spectrum; rather the calculations overshoot the data at 
the nominal p position by about a factor of 2 and under- 
predict the data in the continuum region by about a factor 
of 3 - yet, the integral yield comes out about right. 

In-medium modifications produce a dramatic change: 
both the calculations with an in-medium modified p spec- 
tral function (thick solid line), and with a dropping in- 
medium p mass (thick dashed-dotted line) describe very 
well the marked increase of the continuum yield around 
500 MeV/c^ as well as the depletion at the vacuum p/ui 
position; the differences among the competing approaches 
again are rather subtle. We remark that calculations adopt- 
ing the chiral reduction formalism approach give very sim- 
ilar results, except that the depletion at the vacuum p/oj 
position is absent |115| . 



Before discussing medium modifications in some de- 
tail, let us examine the contributions of p and u to the 
cocktail (for m >400 MeV/c^, say), and more speculative, 
to in-medium pair production. While the uj with its di- 
rect decay and part of its Dalitz decay clearly dominates 
over the p in the cocktail (see Fig. I7.1|l . it is the p which 
provides essentially all of the enhancement. One might 
wonder by which mechanism such drastic change should 
be accomplished. 

Some estimates based on vacuum properties are col- 
lected in the appendix. The yield from mesons formed 
initially by hadronisation is topped by the uj by its larger 
electro-magnetic branching ratio. In contrast, by e^e~ de- 
cays from the in-medium meson population, the p wins 
by a large margin due to its much larger tttt width. Al- 
together, we find that the p/uj ratio of pair yields from 
secondary, in-medium generated mesons, is by orders of 
magnitude larger than for initially produced mesons. 

This digression illustrates from another point of view 
what we know already: the exceptional potential of the 
p propagator to dominate electron-pair production quite 
unrelated to its share in the cocktail - and that there is 
no way to describe the surplus of electron pairs other than 
by direct or thermal radiation out of the fireball, i.e. from 
mesons which are regeneratively produced in tttt annihila- 
tion. The chances are feeble to observe such radiation from 
the Lo, primarily because of its extremely weak coupling 
to the hadronic medium. Medium modifications, however, 
might be observed from primary lo and cf) mesons by the 
fraction that decays within the lifetime of the fireball. 

We return to compare theory to our data. Figure UTTl 
shows that the enhancement at low pair pf^ over the cock- 
tail which is seen in the continuum region and to a lesser 
extent also in the resonance region is a feature present 
in all model calculations, with minor differences only be- 
tween a dropping, broadening, or vacuum p propagator. 
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Fig. 8.2. Comparison of the inclusive mass spectrum of 
Fig. l7.1l to (i) free hadron decays without p decay (thin solid 
line), (ii) model calculations with a vacuum p spectral func- 
tion (thick dashed line) , (iii) with dropping in- medium p mass 
(thick dashed-dotted line), (iv) with a medium- modified p spec- 
tral function (thick solid line). A corresponding cocktail includ- 
ing the p decay is shown in Fig. 17.11 

We meet here the governing features of pion annihilation, 
rather than of medium modifications proper. 

It is the mass spectrum which uncovers the character- 
istic of medium modifications as distinct from pion anni- 
hilation with vacuum p propagator as seen from the com- 
parison of the data with the three modell calculations in 
Fig. 18.21 The differences between the two models incorpo- 
rating medium modifications, however, are rather subtle. 

Let us inspect the model calculations of pf^-selecleA 
mass spectra compared in Fig. 18.31 to the data. The dras- 
tic impact on the shape of the mass spectrum the selection 
of low pf^ has, is also present in the two model calcula- 
tions with modified p. In these models, the effect maybe 
somewhat weaker, yet locally the enhancement over the 
cocktail reaches 10 (see Fig. l7.5l for the complete cocktail). 
The vacuum-p calculation is only weakly affected by the 
P( ^-selection. For larger p'l'^ , data and model calculations 
come much closer to the decay cocktail. 

The processes causing the in-medium changes of the p 
spectral function, or the dropping mass of the p, clearly 
also favour low pair p'^'^ . Such behaviour would arise most 
naturally in the dropping-mass scenario from the Boltz- 
mann (or Bose) factor producing the largest gain for small 
in-medium masses at vanishing 3-momentum (see App. A). 
For the spectral function approach, the observation had 
its impact to install the s-wave N(1520) p-nucleon reso- 
nance as the moving agent in place of the p-wave 
N(1720) resonance which had pioneered the importance 
of p N resonances for medium modifications 117"; this 



change also met requirements by photo absorption data 
to soften the form factor |118j . 

The effective downward shift of strength to lower masses 
in the melting-p treatment |39| is largely due to strong 
meson-baryon coupling, and most approaches agree to its 
importance for generating in- medium effects 27,3 ^115| : a 
finite nucleon chemical potential is required also for some 
meson-meson mixing effects to take place in approaching 
chiral symmetry restoration |119II12U) . Only very small ef- 
fects of baryon density have been reported for UrQMD 
transport calculations ^5,2, •'^"'^ We have remarked that the 
CERES run at reduced SPS energy of 40GeV/n observed 
an even larger enhancement as in 158 GeV/n collisions 
reaffirming the conclusion |121j that the increase in baryon 
density seems to have won over the reduced pion density, 
or lower temperature. As it is the total baryon density 
that matters - vector mesons interact symmetrically with 
baryons and anti-baryons jl22j - the situation at RHIC en- 
ergies will not be greatly different from top SPS energy 
(despite vanishing net baryon density). 

In concluding this review of selected theory descrip- 
tions of the CERES Pb-Au dilepton data, we like to add 
that the spectral function approach had also other suc- 
cesses. Within the same framework, the intermediate-mass 
enhancement observed by NA50 |123| has been success- 
fully described as thermal radiation with a 30% share of 
the quark gluon plasma |124| : there was no need to in- 
voke open-charm enhancement. This finding may be seen 
as the first glimpse of light in the long search for thermal 
qq radiation from the quark-gluon Dlasma [l25lll8) . 

Still within the same framework, the pt spectra of pho- 
tons in 158GeV/n Pb-Pb collisions measured by 
WA98 |12()| have been reproduced; up to transverse mo- 
menta of about 1.5 GeV/c, thermal emission from the 
expanding hadronic fireball has been found to dominate 
with photons mainly of baryonic origin jl27j , fully consis- 
tent with the closely related calculations that describe the 
low-mass dilepton excess observed by CERES. 



9 Conclusion and Outlook 

The analysis of the large unified data sample has sub- 
stantiated the earlier finding of a strong source of con- 
tinuum electron pairs which contributes to the invariant 
mass spectrum beyond the decays of produced mesons, 
most strongly around 500 MeV/c^. There is ample evi- 
dence that we observe dilepton radiation from the interior 
of the hadronic fireball in which pion annihilation medi- 
ated by the p propagator plays a major role. The excess 
yield rises significantly steeper than linearly with charged- 
particle density, consistent with the binary annihilation 
process. Another piece of circumstantial evidence for tttt 
annihilation is delivered by the invariant pair transverse- 
momentum spectra for the continuum pairs of masses be- 

This conclusion rests on the (false) premise that the data 
have been satisfactorily described without medium modifica- 
tions using boosted lo and -q Dalitz decays. 
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Fig. 8.3. Comparison of the pr-selected mass spectra of Fig. [731 < 500 MeV/c (left), > 500 MeV/c (right) to (i) free 
hadron decays without p (thin solid line), (ii) spectral function with vacuum p (thick dashed line), (iii) in- medium dropping p 
mass (thick dashed-dotted line), (iv) in- medium p spectral function (thick solid line). 



tween 200 and 600MeV/c^: the dramatic enhancement 
over the cocktail occurs at very low p^*^ . 

Full calculations with vacuum p describe the measured 
yield about correctly, but fail to account for the charac- 
teristic shape of the spectrum of excess dileptons. The p 
propagator is manifestly modified in the medium which 
is well described by two theories incorporating medium 
modifications of the g, which are, however, very different 
in concept: while the Brown- Rho scaling hypothesis ex- 
plicitly refers to restoration of chiral symmetry, the many- 
body spectral-function approach of Rapp and Wambach, 
although tracing some of the induced mixings, does not 
have chiral symmetry restoration as a ruling concept. That 
it may very well be implicitly included is inferred from 
the fact that the hadronic dilepton rates extrapolated up 
to Tc ('bottom up') come out very similar to the ('top 
down') extrapolated perturbative QGP rates Both 
theories also give a good description of the mass spec- 
tra for selected ranges of pair transverse momentum , 
where an additional preference of the mechanism generat- 
ing the medium modifications has shown up. 

Three unsolved issues remain. To the present data ac- 
curacy, it is not possible to decide between one or the 
other of the two competing theories so that the role of 
chiral symmetry restoration for in-medium modifications 
remains unclear. 

Patience seems also advised in localising the source 
of medium-modified electron-pair production within the 
phase diagram. On one hand, preformation of vector and 
axial-vector correlator strengths in the non-perturbative 
plasma might infiuence dilepton production across the 
phase boundary |128II129| . On the late end of the time 
scale, the impact of a growing pion chemical potential for 
dilepton production might not be entirely settled yet. 



An issue widely overlooked is whether the hadronic 
fireball is 'boiling' long enough to radiate sufficient amounts 
of dileptons. The time spent by the system between chemi- 
cal and kinetic freeze-out came under scrutiny on the basis 
of recent pion interferometric data 11301 which indicated 
only a 30% change in volume. If, in addition, chemical 
freeze-out should occur essentially at the phase boundary 
between hadronic and quark matter, as suggested by the 
asymptotic statistical-model result of Tchem ~ 170 MeV 
at higher energies jlHlj , the purely hadronic origin of the 
low-mass dilepton enhancement might have to be negoti- 
ated again. 

We look out to further progress that can only be ex- 
pected from radically better data, with greatly improved 
statistics and less combinatorial background, but also with 
improved mass resolution. This will not be easy. But a 
first step in this direction is the CERES 2000 run with 
the new Time Projection Chamber; preliminary data have 
been presented very recently |110| . 
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A Appendix p 



decay rate 



The thermal emission rate of dielectrons from two-body de- 
cay of rho mesons has been worked out from Ref. |132| by 
B. Friman and J. Knoll |106| as follows: 



dR 



(l-4m^/M2)3/2 



dMd^q 3(27r)4 (a/2 _ ^2)2 + APr^'^^ 
' M 1 



\/A.f2 + q 2 



(A.l) 



Here, q is the 3-momentum of the rho meson. The new 
expression differs from the one used in the previous GEN- 
ESIS code by the Boltzmann-type phase space factor in the 
curled parentheses being explicitly momentum dependent. To 
obtain the final invariant mass distribution, the differential rate 
d_R/dMd^q has to be integrated over 3-momentum q. The fi- 
nal result used in the 2003 version of New GENESIS [103 for 
comparison to the data reads 



dR a'^rrip 



(1 - 4m2 /m2)3/2 



dM 3(27r)4 (M2 - m2)2 + M^T^^^ 



{2-KMTY'e 



3/2 -M/T 



(A.2) 



B Estimate of p/u ratios 

We estimate the p/uj ratio of electron-pair yields from initially 
produced mesons and from those produced in a hadronic fire- 
ball by using simplified order-of magnitude estimates expressed 
solely by the (vacuum) particle properties |102| . 

Primary mesons are produced during hadronisation. They 
decay into electron pairs with a fraction given by the electro- 
magnetic (e.m.) branching ratio See = ^ee/^- As in p-p col- 
lisions we assume equal primary populations of p and u). 
To the mass range m > 400 MeV/c^, the lo contributes by its 
direct decay and by about 15% of its Dalitz decay. In this mass 
range, the ratio of the number of electron pairs emitted by the 
initial population of p's and tj's is 



Beejp) 



0.3. 



(B.l) 



The p looses in this comparison of e.m. branching ratios due 
its so much larger total width of r{p)= 150 MeV, compared 
to r{uj)= 8.4 MeV. 

We are mostly interested in mesons that decay into electron 
pairs within the lifetime of the fireball since the others can- 
not probe medium modifications. Besides primary produced p 
mesons, an additional source of in-medium decays are from 

Secondary mesons. These are continuously produced in the 
fireball by two-pion annihilation, and most of the time disin- 
tegrate back into two pious as expressed by eqn. (1.1). We like 
to present here a simple estimate of the relative in medium- 
pair production yields from p and uj which neglects all finer 
details of the reaction dynamics by assuming that the ratio 
of rates is given as the ratio of the 2-pion annihilation cross 
sections, which are estimated from detailed balance, times the 
ratio of the e.m. decay widths. The two vector mesons differ in 
an essential manner in strength of coupling to the hadronic 



medium: while the nn channel exhausts the large width of 
the p, -B,r7r(p) ~ 1, the lu (with larger branching into 37r) is 
extremely weakly coupled to the tttt channel, B-mriuj) «2%. 

With Y^ed OC r,r7r -Tee, 



4 X lO"" 



(B.2) 



The estimate demonstrates that by in-medium produced elec- 
tron pairs, the uo meson is completely outnumbered by the p, 
unless its in-medium properties should drastically be changed. 



C Instrumental asymmetries as source of 
enhancement? 

Let us assume there is an asymmetry in reconstruction efficien- 
cies for like-sign pairs (L) as compared to unlike-sign pairs (U) 
which would affect the signal in a most direct way. Asymme- 
tries of this kind might be expected since rings tend to end up 
in a more or less dense environment in RICH-2 depending on 
whether they are deflected away from each other, or not. 

Writing the pair efficiencies as el ~ £o, £u ~ Co (1 + S), 
the signal is given by 



S « ( [/(I - 5) - L)/e = So - 5U/e. 



(C.l) 



To be definite, let us assume a signal So three times the 
hadronic background, i.e. an enhancement factor of three. The 
enhancement will be gone if the asymmetry term reduces the 
signal to one third of So- The corresponding asymmetry pa- 
rameter is 



2(U-L)/3U^-iU- 



L)/L. 



(C.2) 



For a signal-to-background ratio {U ~ L)/L = S/B — 1/13, 
we find that an asymmetry of 5 ~ 5.1% in efficiencies for like- 
sign and unlike-sign pairs is sufficient to fake an enhancement 
factor of 3. 
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